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Abstract 

This paper derives the expressions of the multipole moments of an iso- 
lated gravitating source with an accuracy corresponding to the second post- 
Newtonian (2-PN) approximation of general relativity. The moments are 
obtained by a procedure of matching of the external gravitational field of the 
source to the inner field, and are found to be given by integrals extending 
over the stress-energy distribution of the matter fields and the gravitational 
field. Although this is not manifest on their expressions, the moments have a 
compact support limited to the material source only (they are thus perfectly 
well-defined mathematically). From the multipole moments we deduce the 
expressions of the asymptotic gravitational waveform and associated energy 
generated by the source at the 2-PN approximation. This work, together with 
a forthcoming work devoted to the application to coalescing compact bina- 
ries, will be used in the future observations of gravitational radiation by laser 
interferometric detectors. 
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I. INTRODUCTION 



A. Motivation and relation to other works 

The problem of the generation of gravitational radiation by a slowly moving isolated 
system is presently solved with an accuracy corresponding to the 1.5-post-Newtonian (in 
short 1.5-PN) approximation of general relativity This approximation takes into 

account all contributions in the radiation field at large distances from the system up to the 
relative order e^, where e ~ c^^ is a usual post-Newtonian parameter. The objective of this 
paper is to go one step further in the resolution of this problem by deriving the expression of 
the gravitational radiation field with an accuracy corresponding to the 2-PN approximation, 
taking into account all contributions up to the order e'^. In a forthcoming paper |^ we shall 
consider specifically the case of the radiation generated by a coalescing binary system made 
of compact stars (neutron stars or black holes). 

It is now well-known that post-Newtonian effects in the radiation generated by a coalesc- 
ing compact binary system should be detectable by future gravitational- wave interferometers 
like LIGO and VIRGO [|5HT0|]. It has even been realized in recent years that an extremely 
accurate theoretical signal is required in order to achieve the full potential precision on the 
measurement of the binary's parameters [Q. This remarkable state of affairs is made pos- 
sible by the fact that the signal will spend hundreds to thousands of cycles in the detector 
bandwidth, and that as a result the instantaneous phase of the signal will be amenable to 
a very precise determination. The main influence of the higher-order post-Newtonian cor- 
rections is thus located in the phase of the signal, which is itself mainly determined by the 
rate of decay of the binary's orbit resulting from the gravitational radiation reaction. Note 
that this potential high precision on the extraction of the binary's parameters is interesting 
not only for doing astrophysical measurements of masses, spins, etc. of stars, but also for 
testing some aspects of the nonlinear structure of general relativity ||lT] . 



Analytical and numerical computations of the radiation generated by a small mass in 
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circular orbit around a large one have indicated that post-Newtonian corrections beyond 



the 3-PN approximation may be needed |T^-|T5[. The 2-PN approximation worked out in 



this paper and its sequel thus does not yet reach the ideal precision required by the 
observations of coalescing binaries, but still represents an appreciable improvement of the 
existing situation. 

The investigation of this paper will be based on a "multipolar-post-Minkowskian" for- 
malism for dealing with the gravitational field in the external vacuum region of the source. 
Such a formalism combines a multipole moment expansion, valid in the exterior of the source, 
with a post-Minkowskian expansion, i.e. a nonlinearity expansion or expansion in powers of 
Newton's constant G, valid wherever the field is weak. This double expansion formalism is 
originally due to Bonnor and collaborators ||16|-p!9[], and has been later refined and clarified 



by Thorne ||20|. In recent years, the formalism has been implemented in an explicit and con- 



structive way pT] , p2| (with the help of some mathematical tools like analytic continuation), 
and applications have been made to the problems of gravitational radiation reaction p3|j2^ 
and wave generation . The field being determined only in the exterior of the source, the 
multipolar-post-Minkowskian formalism must be supplemented by a method of matching to 
the field inside the source. We shall use a variant of the well-known method of matching of 
asymptotic expansions (see e.g. |^) which has already, on several occasions, found its way 
in general relativity [^^-|29|. 

Inspection of the solution of the wave generation problem at the 1.5-PN approximation 
readily shows what is needed for extension at the 2-PN approximation. Indeed, the 
works in and have obtained respectively the expressions of the mass-type and current- 
type multipole moments of the source with relative precision e^, the next-order correction 
being of order 0{£^). The inclusion done in |I| of the terms in the mass multipole moments 
permits to solve the wave generation problem at the 1-PN approximation (in fact, only the 
terms in the quadrupole mass moment are needed). This has set on a solid (and well- 
defined) footing previous works by Epstein and Wagoner [Q, and Thorne who obtained 



formally correct but divergent expressions of the moments at the 1-PN approximation (see 



|T|] for discussion). In the work the nonhnear effects in the radiation field were added 
and shown to arise at the level e^, essentially due to the "tails" of gravitational waves. 
Higher-order nonlinear effects were shown to be of order 0{e^) at least. Now recall that the 
contribution of a moment with multipolarity / scales like in the radiation field, and that the 
current moments always carry an additional factor e as compared with the corresponding 
mass moments. Therefore we conclude that what is needed for solving the 2-PN wave 
generation problem is only to find the extension of the expression of the mass quadrupole 
moment of the source with relative precision e^. Note that the computation of the moments 
of the source can be done in the near zone of the source. 

This paper will thus mainly focus on the matching between the external and internal 
gravitational fields in the near zone of the source (both fields are expressed in the form of 
a post-Newtonian expansion when £ — > 0) with a precision consistent with the inclusion of 
the terms in the mass multipole moments. (Although the quadrupole mass moment is 
sufficient for our purpose, we shall compute the terms in all mass moments of arbitrary 
multipolarity /.) The matching performed in [|T| was based on a particularly simple closed 
form of the internal gravitational field of the source including 1-PN corrections. On the 
other hand, the matching performed in made use of some specific distributional kernels 
to deal with the quadratic nonlinearities of Einstein's equations. None of these methods 
can be applied to our problem, which necessitates the inclusion of the full 2-PN corrections 
in the field, depending not only on quadratic but also on cubic nonlinearities of Einstein's 
equations. In this paper we shall employ a matching procedure which is far more general 
than the ones followed in . In particular, we shall show how one can deal with the cubic 
nonlinearities of Einstein's equations without any use of distributional kernels. 

The end result obtained in this paper expresses the multipole moments of the source as 
integrals extending over the distribution of stress-energy of the matter fields in the source 
and of the gravitational field. This result is similar to the one we would obtain by using as 
the source of the radiation field the total stress-energy (pseudo) tensor of the matter and 
gravitational fields, and by considering formally that this tensor has a spatially compact 
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support limited to the material source. It is well-known that by proceeding in this formal 
way (i.e., d la Epstein- Wagoner-Thorne), we obtain integral expressions of the multipole 
moments which are divergent, because of the non-compact-supported distribution of the 
gravitational field. (Indeed, the integral expression of a moment with multipolarity / involves 
a large power ~ |x|' of the radial distance to the source, which blows up at infinity.) We prove 
in this paper that the correct expressions of the multipole moments must involve also an 
analytic continuation factor |x|^, where S is a complex number, and a "finite part at B = 0" 
prescription to deal with the a priori bad behavior of the integrals at spatial infinity. In 
such a way the expressions of the moments are perfectly well-defined mathematically. When 
a pole ~ 1/5 occurs in an integral due to the bad behavior of the integral at spatial infinity, 
the "finite part at B = 0" introduces an additional contribution to the moment which must 
absolutely be taken into account. (However, we shall see that no poles occur at the level 
investigated in this paper.) 

Let us stress that the expressions presented here of the multipole moments are not 
manifestly of compact-supported form (contrarily to say the lowest-order expressions at 
the 1-PN level [Q), but are numerically equal, thanks to the properties of the analytic 
continuation, to some compact-supported expressions which could be constructed at the 
price of introducing more complicated potentials besides the usual Poisson integrals of the 
mass and current densities in the source. Such a construction is however unnecessary (and 
is somewhat awkward) in practical computations of the moments for specific sources. 

It is likely that the expressions of the moments expressed in this way as integrals ex- 
tending over the total stress-energy distribution of the matter and gravitational fields and 
regularized by means of analytic continuation, will admit a generalization to higher non- 
linearities in the field, and higher post-Newtonian approximations. Such a generalization, 
which will be the subject of future work, should permit the resolution, at least in principle, 
of the problem of the generation of gravitational radiation by a general isolated system up to 
the high level of approximation required by the observations of coalescing compact binaries. 

This paper is organized as follows. In section |l| we compute the gravitational field 
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both in the near zone of the source (where a direct post-Newtonian expansion of the field 
equations is performed), and in the external near zone of the source (where we use the 
multipolar-post-Minkowskian solution of the vacuum equations). In section ^ we impose 
that the two fields are isometric in their common domain of validity, namely the external 
near zone of the source. This yields a matching equation valid up to an appropriate post- 
Newtonian order, and which is used in section to obtain the expressions of both the mass 
and current multipole moments of the source. The formulas for the asymptotic waveform 
and the associated energy generated by the source at the 2-PN approximation are also given 



in section IV. The paper ends with three appendices, and we start with our notation for 



Einstein's equations. 

B. Notation for Einstein's equations 

Throughout most of this paper, we shall use Einstein's equations in harmonic coordinates. 
The field deviation from Minkowski's metric is denoted by /i"^ = yZ—gg^^ — f]"^, where greek 
indices a,l3,fi, v... range from to 3, gap is the usual covariant metric, g""^ is the inverse of 
ga0, g is the determinant of gap-, and 77°^ = diag{—l, 1, 1, 1) is the Minkowski metric. The 
condition of harmonic coordinates reads as 

dph''^ = , (1.1) 

where dp = d/dx^ is the usual derivation with respect to the harmonic coordinates. Ein- 



stein's equations reduced by the condition (|1.1| ) are written in the form 



□/,a/3^1MAr-'3 + A"/3(/,) . (1.2) 



We denote by □ = rj'^^dadp the flat d'Alembertian operator, by the stress-energy tensor 
of the non-gravitational fields (T"^ has the dimension of an energy density), and by A the 
absolute value of g, i.e. X = \g\ = —g. In terms of a series expansion in the field variable 
fi^is have 
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X = l + h + -ih^-h^,h^'') + 0{h^) , 



;i.3) 



where h^^^, = 'q^aVup^"^ h = "qaph^^ . The second term in ( p..2|) is an effective gravita- 
tional nonhnear source including all the nonlinearities (quadratic, cubic, . . . ) of Einstein's 
equations. We denote 

K^P[h) = N'^'^ih, h) + M"'3(/i, h, h) + 0(/i^) 



:i.4) 



where the quadratically nonlinear term is given by 



+r] 



a(3 



4 o Z _ 



;i.5) 



and where the cubically nonlinear term is given by 

M"^(/i, h, h) = -h^'^i^d^h^pd^hP^ + dphyphl - dphy^hf^p) 



(^d^^Kpd^h-p - 2d,h^p>d^h''p - U^^hd^h^ 



al3 



h^'dphd^h - h'^'^dphp^dPh - ^hP'^dphpM'"' 
1 



^hP'^d^Kd'^hP + -hP'^d.KdPK 



In (|1.5|) and ( p..6|) , we raise and lower all indices with the Minkowski metric, and we denote 
t{ai3) = + tpo). By taking the divergence of ( |1.5| ) and (|1.6|), we obtain the relations 

dpN^^ = [-dpK + V - \d''h r/^.) W , (1.7) 



+ 



-h'^Pdphp, - hppd'^K + K.d^h'^p 



+\hp,d''h + hp^d'^hP'^r^p, - h"Pdphr]p, 



ahp" 



which are consistent, by Bianchi's identities, with the conservation (in the covariant sense) 
of the stress-energy tensor of the matter fields. 
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II. THE INTERNAL AND EXTERNAL GRAVITATIONAL FIELDS 



A. Solution of Einstein's equations in the internal near zone 

Let us define the internal near zone of tlie source to be a domain Di = {(x, t)/|x| < r^} 
wliose radius rj is adjusted so tliat (i) rj > a, wliere a is tlie radius of a sphere which totally 
encloses the source, and (ii) rj -C X, where X ~ ac/v is the (reduced) wavelength of the 
emitted gravitational radiation, and f is a typical internal velocity in the source. 

Defining Di in such a way assumes in particular that a <^ A, or equivalently £ ^ 1 
where e ~ f/c is a small "post-Newtonian" parameter appropriate to the description of 
slowly moving sources. (We shall also assume that the source is self-gravitating so that 
GM/ (ac^) ~ e"^ where M is the total mass of the source, and that the internal stresses are 
such that T^i jT^^ ~ e^.) In Dj, we can solve Einstein's equations ( |l.l|) -( pr2|) by formally 
taking the post-Newtonian limit e ^ 0. The following notation will be used for terms of 
small order in the post-Newtonian parameter e. By A = 0{p) we mean that A is of order 
0{e'^)\ by 5" = 0{p,q) we mean that the zero component of the vector B°' is = 0{p), 
and that the spatial components of B°' are B^ = 0{q); and similarly by C"^ = 0{p, q, r) we 
mean that C^^ = 0{p), = 0{q) and C*-' = 0(r). (The latin indices . .range from 1 
to 3.) 

From the contravariant components of the stress energy tensor of the matter fields, 
we define a mass density a, a current density cTj and a stress density aij by the formulas 

^00 _j_ ^ii 

a= , (2.1a) 

a. = — , (2.1b) 

c 

(y^,=T^ , (2.1c) 



where in (|2.1a| ) = 115ijT^^ denotes the spatial trace of T"^. In all this paper, we shall 



assume that the matter densities cr, cxj and aij are of order e° when e — >• 0, i.e. 



a, (Ji, (Jij = 0(0) . (2.2) 



We introduce next some reiardec? potentials generated by the densities a, ai and aij. First, V 
and Vi are the usual retarded scalar and vector potentials of the mass and current densities 
0" and CTj, i.e. 

y(x, t) = G [ r^-^ a(x', t - i|x - x'l) , (2.3a) 
J |x — x'l c 

y,(x,t) = g| a,(x',t- i|x-x'|) , (2.3b) 

satisfying nV = —AnGa and = —AnGai. Second, Wij is a more complicated retarded 
tensor potential defined by 



iy,,(x,t) = G j 



X — X' 



(x',t--|x-x'|) . (2.4) 



c 



Note that while the potentials V and have a compact support limited to that portion of the 
past null cone issued from the field point (x', t) which intersects the source, the potential Wij 
has not a compact support. From ( |2.3a| ), we see that r(x',i(;-i|x-x'|) behaves like GM/|x'| 
when |x'| oo, where M = J d^yaly, — oo) +0(2) is the initial mass (or ADM mass) of the 
source, and we can check from this that the potential Wij is given as a convergent integral. 
We shall often abbreviate formulas such as ( p.3| - p.4D by denoting the retarded integral of 
some source /(x, t), having adequate fall-off properties along past null cones, by 

(□^V)(x,t) = --L / ^!!Z^/(x',t-l|x-x'|) , (2.5) 
An J |x — x'l c 

so that e.g. V = -AnGa~^a. 

It is convenient in most of this paper to keep the potentials (p.3[ - |2.4| ) in retarded form, i.e. 



to not expand when c — * +oo the retardation argument t — |x — x'|/c. This permits to avoid 
possible problems of convergence with the expansion of potentials with non-compact support. 
However, most of the equations of this paper will be only valid up to some remainder in 
the expansion e ~ c~^ — 0, and we can replace, when there is no problem of convergence, 
the retarded potentials by their post-Newtonian expanded forms up to the accuracy of the 
remainder. For instance, from the Newtonian equations of continuity and of motion. 
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dta + diai = 0{2) , (2.6a) 
dtai + djai, = ad^V + 0(2) , (2.6b) 

we deduce that the potentials V , Vi and Wij satisfy the conservation laws 

dtV + d,V, = 0{2) , (2.7a) 
dtV, + d,W,, = 0{2) . (2.7b) 

In (|2.6b|) and (|2.7| ) we can replace the potentials V , Vi and Wij by corresponding Poisson- 
type potentials, e.g., V = U + 0(2) where U is the Newtonian potential of the mass density 
cr, satisfying Af/ = —AiiGa. 

We now proceed to solve Einstein's equations ( [1.1H1.2| ) with an accuracy corresponding 
to the post-Newtonian order 0(8, 7, 8), by which we mean 0(8) in the 00 and ij components 
of /i"'^, and 0(7) in its Oi components. The insertion of the lowest-order results h^^ = 
-4\//c2 + 0(2), h^' = 0(3) and h'^ = 0(4) into the right-hand-side of ([l^) with the exphcit 
expression ( |1.5| ), yields first the equations to be solved at order 0(6, 5, 6). We get 

Oh'' = ^ (l + 5") - ^^'^^^'^^ + ^(6) , (2.8a) 
ah'^ = l^T'i + 0{5) , (2.8b) 
ah'J = ^r^' + ^ [d,Vd,V - ^6,,dkVd,V^ + 0(6) . (2.8c) 

These equations can be straightforwardly solved by means of the potentials V, Vi and Wij 
we defined in (|2.3|) and ( ^^ , and by means of the trace W = Wa = T,6ijWij of the potential 



Wij. The result is 



= + ^{W- 2V') + 0(6) , (2.9a) 

h^^ = -^V, + 0{b) , (2.9b) 
h^j = -^Wi, + 0(6) . (2.9c) 



Since the potentials satisfy the conservation laws ( |2.7| ), we see that the field ( |2.9| ) satisfies 
the approximate harmonic gauge condition dph'^^ = 0(5,6). 
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The next step consists in iterating Einstein's equations by using the expression (|2.9|) of 
the field. We substitute ( p.9| ) into the right-hand-side of ( |1.2| ) with the help of ( |1.3| ) and 
of the explicit expressions ( |1.5| ) and ( |1.6| ) of the quadratic and cubic nonlinearities. In this 
way, we find the equations to the satisfied by the field up to the order 0(8, 7, 8), which are 

Uh^P = ^^KV, W)T''f + T\V, V„ W,,) + 0(8, 7, 8) . (2.10) 

In these equations, we denote by A the post-Newtonian expansion of X = \g\ = —g when 
truncated at the order 0(6), i.e. A = A + 0(6), and expressed with the potentials V and 
W = Wii. From (|OI) we have 

\{y,W) = l + \v --AW -V^) . (2.11) 

(Note that because the matter stress-energy tensor is T"^ = 0(-2, -1, 0), the precision on A 
as given by ( p. 11 ) is necessary only in the 00 component of the equations (|2.10|) .) Similarly, 



A*^^ in ( 2.101) is defined to be the post-Newtonian expansion of the effective nonlinear source 



in the right-hand-side of Einstein's equations (|1.2| ) when truncated at order 0(8,7,8), i.e. 
A"^ = A"^ + 0(8, 7, 8), and expressed in terms of the potentials V , Vi and Wij. From ( p..5|) 
and (|1.6| ), we arrive at the expressions 

-Wkmdty + \{dtVf + \duV^{dkVm + 39^\4) 

o Z 



+dkVdtVk + 2dkVdkW - -vdkVduV^ , (2.12a) 

A°*(y, F„ W,j) = — \dkV{d^^k - dkV,) + -dtVdiV^ , (2.12b) 
V\V, V,, W,j) = ^ S^d,Vd,V - h,,dkVdkV^ + ^l^2d^,VdtV,) - d,Vkd,Vk 

-dkVAV, + 2d(^,VAVj) - h,{dtVf - 6,,dkVdtVk 

o 

+ls,,dkVUdkVm - dmVk)^ . (2.12c) 

By using the equations of motion and of continuity including e"^ terms (which are given by 
( pSD in the appendix B below), one can check that the expressions (|2.11| ) and (p.l2|) imply 
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0(7,8) . 



(2.13) 



Note that the condition ( p.l3|) checks all terms in ( p.l2|) except the terms in the 00 
component of (|2.12a|) and the term in A. We have checked these terms by adding 
in the Oi component A*^* of (|2.12b|) the next order c"'' terms and computing the divergence. 
Note also that the consideration of the next order c terms in A would allow the control 
of the correction terms not only in the mass-type source moments (which is the aim of 
this paper), but also in the current-type source moments. However, we have chosen here not 
to include these terms because they are not necessary for solving the 2-PN wave generation 
problem, and because they somewhat complicate the discussion with the need of introducing 
new potentials besides V, Vi and VFjj, and the necessity of a better control of the external 
metric in § |IIB| . These terms will be considered in a future work where we investigate how the 
procedure followed in this paper could be systematically extended to higher post-Newtonian 
orders. 

Finally, a solution of ( p.lOj) , with the required precision, can simply be written as 

h-P = [1^AT°^ + r^l + 0(8, 7, 8) , (2.14) 



where is the retarded integral operator defined in (|2.5| ). By ( |2.13| ), this solution satisfies 
also the approximate harmonic gauge condition 

= 0(7,8) . (2.15) 

Note that we could have added in ( |2.14| ) some homogeneous solutions of the wave equation 
which are regular in Di and satisfy the harmonic gauge condition ( |2.15| ). It is simpler to 
use the solution ( p.l4| ) as it stands since we shall show that it directly matches the exterior 
metric. 



B. Solution of Einstein's equations in the external near zone 

Let = {(x, t) / r = |x| > Te} be an external domain surrounding the source, where 
Te is adjusted so that a < < rj, a being the radius of the source and rj the radius of 

13 



the inner domain Di defined in §11 A| . By our assumption GM/{ac^) ~ e^, gravity is weak 



everywhere and in particular in Dg, so we can solve Einstein's vacuum equations in D^. 
by means of the multipolar and post-Minkowskian approximation method developed in our 
previous works pT|-p^ on foundations laid by Bonnor ||I6| and Thorne . More specifically, 



we use the construction of the external field which is defined in § 4.3 of |2T|] and which is 
referred there to as the "canonical" external field. Using a formal infinite post-Minkowskian 
expansion, or expansion in powers of Newton's parameter G, the "canonical" external field 



KL = V-fcan g^^n " IS glVeU aS 

= G Km,^ + G'Km,^ + ... + G-hZin) + ■■■, (2.16) 

where the coefficients h'^'^^^^^) arbitrary nth power of G are algorithmically constructed 

from the knowledge of the previous coefficients h^f^^^-^ (with m < ri) by post-Minkowskian 
iteration of the vacuum equations. This iteration is made possible by a systematic use of 
multipole expansions, valid outside the source (in De), for all the coefficients /ican(n) ( P-16D . 



The whole iteration in ( |2.16D rests on the first of the coefficients, namely the "linearized" 



field /ican(i), which is chosen to be the most general solution, modulo an arbitrary linear 
gauge transformation, of the linear vacuum equations. This solution reads as 

G0) = -^^'^* ' (2-17a) 
G/^L(i) = -4v;^^\ (2.17b) 
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where the external potentials V^'^^*, V^^^^ V^^^* (different from the inner potentials V , Vi, 
Wij) are given by some explicit infinite multipole expansions of retarded spherical waves 
(solutions of the homogeneous wave equation in Dg), namely 

"1 - . / r' 



t^o ^! Lr V cj 



(2.18a) 



iL-l , 

r V c 
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aL~l 



s, 



r 



bL-l 



cxt 



1=2 



-d 



L-2 



(1) 

j)bL~2 



r 

t - - 



(2.18b) 



(2.18c) 



See e.g. (8.12) in Thorne pO|. Our notation is as follows (anticipating also on future needs). 
Upper case latin letters denote multi-indices with the corresponding lower case letters being 
the number of indices, e.g. L = iii2 ■ ■ - ie- Similarly, L — l = ii---i£_i,L — 2 = ii--- i^_2 and 
aL = ail ■ ■ -ii. A product of space derivatives di = d/dx^ is denoted by di = di^di^ ' ' '^ir 
Similarly, = ■ ■ ■ x^'^ and = n^^rt'^ ■ ■ ■ n** where n* = xY|x| = x^ /r. The symmetric 

and trace-free (STF) projection is denoted with a hat, e.g. di or x^, or sometimes by e.g. 
di^L)- M^P\t) denotes the pth time derivative of M{t), and Tij = |(Tjj + Tji). 



The linearized external field ( p.lTHzT^ ) satisfies the linearized equations ^h'^!^^ 



'can(l) 







everywhere except at the spatial origin r = of the coordinates. The harmonic gauge 



condition dyh'^^^(^-^ = follows from the (exact) identities 



dtVr' + o'.V^.f * = . 



(2.19a) 
(2.19b) 



Note that the potential Vf^^ is trace-free: 



ext 







(2.20) 



As we see, the potentials ( [^.18| ) depend on two infinite sets of functions of time, Miit) for C = 
0, • • ■ , cx), and Siit) for £ = 1, ■ ■ ■ , oo. These functions are STF in their i indices. They can 
be viewed respectively as some "canonical" mass-type and current-type multipole moments 
parametrizing the external canonical metric. They are completely arbitrary functions of 
time except that the lowest multipole moments M (mass monopole). Mi (mass dipole) and 



Si (current dipole) are constant: M^^^ = M^^^ = S^^^ = 0. Note that it was assumed in 



:(i) 



2l| that the moments ML{t) and Siit) are constant before some remote date in the past. 
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We shall admit here that one can cover a more general situation where (for instance) the 
ith time-derivatives of M^lt) and 5'i(t) become asymptotically constant when t — oo, 
corresponding to a situation of initial scattering in the infinite past. Furthermore, we shall 



relax without justification the assumption made in [^] that the multipole expansions are 
finite; that is, we assume that we really have two infinite sets of moments M^lt) and Silt). 
These two amendments almost certainly have no incidence on the results derived in this 
paper. 

Starting with the linearized metric ( p.l71 - |2!T^) , one constructs iteratively the nth coeffi- 



cient /?-can(n) series (|2.16|) by the formula 

Cn(n) = FPb=0 [r^'AZin)] + Cn(n) • (2-21) 

The first term in this formula involves A'^!' / ^ which is defined to be the coefficient of G" in 
the expansion of the effective gravitational source A.^^{h) in the right-hand-side of Einstein's 
equation (|1.2| ) and computed with the "canonical" field (|2.16|) . That is, 

A^'^l/^can) = G'KLm + ■ ■ ■ + G-A^^r.H + ■ ■ ■ . (2.22) 

Since is at least quadratic in Aclnin) ^'-'^ n is a function only of the previous n — 1 
coefficients /ican{i)5 ■ ■ ■ , ^can{n-i)- For instance, we have (with evident notation) 

A^:n(2)=^'^'(W)>W))> (2-23a) 

^^11(3) = A^'^''(^can(l), ^can(2)) + A^'''' (^can(2) , ^can{l) ) + ^^"^ (^can(l) , ^can{l) , ^can(l) ) , (2.23b) 



where N'^^ and M^'^ are given by ( |1.4|) -( pr6|) . The retarded integral operator D/j^, which is 



defined in ( p.5| ), acts on the source A^^^^^^^ but multiplied by an analytic continuation factor 
r^, where r = |x| and i? is a complex number. The introduction of this factor is required 
because the linearized metric /ican(i) of (|2.17| )-( |2rT8| ), and all subsequent metrics /ican(n) and 



sources Acan(n); circ valid only in the exterior of the source (in D^) and are singular at the 
spatial origin of the coordinates, r = 0, located within the source. It has been shown in 
||2T| that for B a complex number the retarded integral f{B) = □]j^(r^Acan(n)) defines an 
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analytic function of B all over the complex plane except in general at integer values of B. 
Near the value B = 0, f{B) admits a Laurent expansion of the type f{B) = HapB^, where 
p E The coefficient of B^ in this expansion, i.e. oq, is what we call the finite part at 
B = (or FPb=o) of the retarded integral f{B). This is the first term in ( p.21[ ); it satisfies 
□oq = Acan(n) (and is also singular at the origin). Thus the introduction of the analytic 
continuation factor is a mean (and a convenient one) to obtain a solution of the Einstein 
equation (|1.2| ) in — this is the only thing we need (see for more details about this way 
of proceeding). The second term gcan(n) in (|2.21| ) is a particular retarded solution of the wave 



equation, i.e. □gcan{n) = 0, whose divergence is the opposite of the divergence of the first 
term, and thus which permits to ensure the satisfaction of the harmonic coordinate condition 
^''^can(n) ~ 0. The prccisc definition of the term gcan(n) is reported in the appendix A where 
we control its order of magnitude in the post-Newtonian expansion. 

The external field ( |2.16|) , in which we have (|2.17|) - (|2.18|) and (|2.21|) , represents the most 



general solution of Einstein's equations in De, and is parametrized by the arbitrary "canoni- 
cal" multipole moments ML{t) and S'i(t). Now the point is that one knows (see P^P],|3|) how 
to relate the observable "radiative" multipole moments UL{t) and ^^(t), parametrizing the 
field in the distant wave zone of the source (where the detector is located), to the moments 
Miit) and Siit) (see § |rv| below). Therefore, what is only needed in order to compute the 
distant wave field is to consider the near-zone expansion, or expansion when c +00 or 
e — s> 0, of the external field (|2.16|) up to some suitable order, so as to give by matching to 
the inner metric constructed in § II A a suitably accurate physical meaning to the moments 
Mi(t) and SL{t) in terms of the source's parameters. Equivalently, this means that we must 
consider the external field ( p.l6|) in that part of the external domain which belongs to 
the near zone, i.e. D^n Di = {(x, t)/re < r < r^}, in which we can simultaneously expand 
the external field when e and keep its multipole moment structure. First of all, from 
the fact that an arbitrary nonlinear coefficient /ican(n) ( |2.16| ) is of order 0{2n, 2n + 1, 2n) 



when 6^0 (see (5.5) in |]2T]]), we find that the neglect of all nonlinear iterations with n > 4 
permits the computation of the field up to the order 0(8, 9, 8) when e ^ 0, i.e. 
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K^n = GKmi) + G^Cn(2) + G^CnO) + ^(8, 9, 8) . (2.24) 
Secondly, it is shown in the appendix A that the second terms g^an(2) '?caii(3) 



definitions of the quadratic and cubic coefficients (|2.21|) have (at least) the following orders 
of magnitude when e — > 0: 

Cn(2) = 0(7,7,7), (2.25a) 
Cn(3) = 0(8,7,8). (2.25b) 



Therefore, from (|2.21|) , we can write the expansion when e — of the canonical external 



field ( p.24| ) in the form 



hi^L = GKU) + FFb=o O-^' \r''{G'AZi2) + G'AZis))] + 0(7, 7, 7) , (2.26) 



(2) ^ "'-can(3)/ 



where the remainder term 0(7, 7, 7) is dominated by the contribution ( p. 25a)) of q'^an{2)- 



The source terms Acan(2) and Acan(3) in (|2.26|) are given more explicitly by ( |2.23|) . To 



compute their sum with an accuracy consistent with the remainder in ( 2.26 ), we must know 
the quadratic metric /ican(2) up to the order 0(6,5,6). This is of course analogous to the 
computation we have done in § IIA, where we had first to solve Einstein's equations up to the 
order 0(6, 5, 6) before reaching the looked-for order 0(8, 7, 8). The quadratic source at the 
order 0(6,5,6) is readily obtained by substituting (|2.17|) - (|2.18|) into (|2.23a|) and discarding 
0(6,5,6) terms. We obtain 



,^2 A 00 

^^can{2) — 


_llafci^^=^*afc\/^^* + o(6) , 


(2.27a) 


^2 A Oi 

^^can{2) — 


0(5) , 


(2.27b) 


/^2 A ij _ 
^ ^^can(2) - 


^ l^d.V'^'d.V'''' - ^S.jdkV'^'dkV'''''^ + 0(6) , 


(2.27c) 


we deduce 








GXM2) = -^AV^^r + Oi6), 


(2.28a) 




G'^h^ln{2) =0(5) , 


(2.28b) 




0^/^1(2) = -^^? + 0(6), 


(2.28c) 



where we have introduced the new external potential 



, . ............. , 1 



Zlf = FPR=n a-^ 



(2.29) 



The justification of ( p.28| ) is as follows. We know that the post-Newtonian expansion of 
the regularized retarded operator FPD^^ acting on terms belonging to the quadratic source 
^can(2) is equal to the expansion obtained from the action of the regularized "instantaneous" 
operator FP/"^ := FPS^o(^/c5^)^''^"^"S where A'^^-^ is the (A; + l)th iteration of the 
Poisson operator A^^, modulo negligible terms of order 0(10, 9, 8) (see e.g. (3.7) and (3.21), 
with n = 2, in [^). Thus, FPD^^ acting on the remainder 0(6, 5, 6) in ( p.27|) is itself of order 



0(6, 5, 6). Furthermore, we know that FPA ^ (namely the first term in FPJ ^) gives simply 
j^^cxt^2 T^\yQYy acting on A[([/''^*)^], where f/'^^* is the "Newtonian" potential associated with 
and given by f/^^* = GS^((9Lr-i)Mi(t) (see e.g p. 395). Since l/'^^* = f/^^* + 0(2), 
we have FPn;^^(A[(\/'=^*)2]) = (F^^t)2 + 0(2) as has been used in ( p:28a| ). Note that this 



last fact implies that the trace Z''^'^ = Z^f" of the potential ( f^.29[ ) satisfies 



= l(\/^-t)2 + 0(2) . (2.30) 



Using now ( p.30| ) and ( |2.20|) , we can write the canonical field up to the order 0(6, 5, 6) in a 



form which is formally identical to that of the inner field (|2.9| ), i.e. 



/^oo^ = -1 + 1 (1^-t - 2(V'=^*)2) + 0(6) , (2.31a) 
^cL = -4^r* + 0(5), (2.31b) 



4 



hln = --,W^f + 0{Q) , (2.31c) 



where we have defined 



W^f = V^,f * + Zif , (2.32a) 



whose trace is 



p^ext ^ ^ext _ ^cxt _ 1(^^x1^)2 ^ ^^2) . (2.32b) 
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Finally, since the "canonical" field ( p.31| ) has the same form in terms of the external 
potentials as the inner field ( |2.9| ) has in terms of the inner potentials, it is clear that the 
source term G'^K'^l^^^) + ^'^^crn(3) right-hand-side of ( p.26| ) will be equal modulo 

the same order 0(8,7,8) as in ( p.lO| ) to the truncated source defined in (^), but 
computed with the external potentials V*^^^ , V^^^ and W^^^^ instead of the inner potentials 
V , Vi, Wij. (Indeed, never a Laplace or d'Alembertian operator enters in the effective 
gravitational source of Einstein's equations, which would give a different result when acting 
on an external potential or on an inner one.) Thus, we can write 

G'Km2) + G'Km3) = A'\V^^\ Vr\ W:f) + 0(8, 7, 8) , (2.33) 

where the right-hand-side is obtained from ( p.l2|) by the simple replacement V, Vi, Wij 
^ext^-^ext^™t_ Inserting (p33|) into and recalling that FPD^^ acting on 0(8,7,8) 



is also 0(8, 7, 8), yields our looked- for expression of the external field, namely 



h^L = GhZw + FPb=o r^A^'^F-*, Vr\ + 0(7, 7, 7) . (2.34) 



1(1) 



B \ 1^^^ / T rcxt T^ext TI/ext^ 



It is to be noticed that the remainder in ( |2.34D is 0(7, 7, 7) instead of the remainder 0(8, 7, 8) 



in the inner field ( p^.l4D because it involves the not controlled contribution 0(7, 7, 7) of g, 



flU 

can(2) 



in ( |2.25a|) . In the next section we match the external field ( p.34|) to the corresponding inner 



field (|1D 



III. MATCHING OF THE INTERNAL AND EXTERNAL FIELDS 



A. Coordinate transformation between the internal and external fields 



We require that the internal field h'^^ constructed in § [II A| and the external field hf^!^^ 



constructed in § [II B| are isometric in their common domain of validity, which is the exterior 
near-zone domain Di (1 De = {(x, )f:)/re < r < Vi} of the source. Denoting by x" the 
harmonic coordinates used in the inner domain Di (see § [II A[), and by x^^^^ the "canonical" 
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harmonic coordinates used in the exterior domain (§ |IIB| ), we thus look for a compatible 
coordinate transformation 



x^,,^{x)=x^ + ^^{x) , (3.1) 
where the vector (f'^{x) is assumed to be in the form of a multipolar and post-Newtonian 



expansion appropriate in Di f] De- [In § |IIB| , we have for notational convenience abusively 



denoted by x'^ what really are the canonical harmonic coordinates x^i^^.j Since the two 
coordinate systems x'^ and are harmonic, the vector ip^ satisfies the (exact) relation 

+ h-^{x)dl^^^ = , (3.2) 

where □ = fj'^^d'^p. The (also exact) transformation law of the field deviation /i"^ under the 
change of coordinates ( p.l|) is given by 

V'-" + KLi^c^r.) = ^(5^^ + d^^niS; + dp^niv"-^ + h-^ix)] , (3.3) 

where J = det(9a;can/c^a^) denotes the Jacobian determinant of the coordinate transforma- 
tion. 

We now expand the transformation law ( p. 3] ) when £ — > 0. The field deviation h°'^ in 
the right-hand-side of ( p.3|) is by ( p.9|) of order 0(2,3,4) when e 0. Furthermore, let us 
assume that the vector yj^ in the coordinate transformation (|3.1|) is of order 



^^ = 0(3,4). (3.4) 

This assumption (which has already been made in previous works [0,0]) is proved below 
when we show that it leads to a consistent matching. We can then easily see that the 
transformation law ( |3.3| ) reduces up to the order 0(6, 7, 8) to a linear transformation, 

K^ni^) = h^\x) + d^^^'ix) + 0(6, 7, 8) , (3.5) 

where we have expressed both sides of the equation in terms of the inner coordinates x^, 
and where dip^^ denotes the linear part of the coordinate transformation given by 
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diff''' = d^'ip" + - r]''''dxip^ . (3.6) 

By taking the divergence of ( p.6|) and using (|3.2|) we obtain 

d^d^^-" = = -hP'^dl^if^ , (3.7a) 

from which we deduce the order of magnitude 

d^d^"" = 0{7,8) . (3.7b) 

Note that one could have a priori expected some nonhnear terms to appear at the order 
in the ij component of the equation ( |3.5| ) — for instance, terms hke h'^^d^^ip^^ or h^^'^d^^Lp^. 
Such nonhnear terms are, however, absent at this order (as was also found in 0). On the 
contrary, some nonlinear terms arise at the order in the 00 component of (|3.5|). These 



terms, which will be needed in the following, are obtained by a short computation showing 
that the 00 component of ( |3.5|) , now valid up to the order 0(8), is given by 



/iOO^(a;) = /lOC(x) + V° + 2/i°^9^v^° - 9^(/i°V) + 5i^°5i<^° + 0(8) . (3.8) 

Here also we have expressed both sides of the equation in terms of the inner coordinate 
system x^. 

B. Matching of the compact-supported potentials V and Vi 

A requisite in the matching procedure is to find the relations linking the multipole ex- 
pansions outside the source of the inner potentials V, Vi and Wij and the external potentials 
p^ext^ l^ext g^j-^^ Wij^^. We deal in this subsection with the case of the compact-supported 
potentials V and Vi. The more complicated case of the non-compact-supported potential 
Wij is reported in the next subsection. 

The inner and outer fields have been shown in § |I|to take, up to the order 0(6, 5, 6), the 
same functional forms ( p^.9D and ( |2.31| ) in terms of their respective potentials. On the other 



hand, the coordinate transformation between them takes, up to this order (and even up to 
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a higher order), the hnear form ( |3.5| ) with ( |3.6| ). We first substitute into the sum of the 00 
component and of the spatial trace ii of (|3.5| ) the inner and outer fields ( |2.9| ) and ( ^.311 ). 
This leads to an equation whose solution is easily seen to be 



^ext = v + cdtip^ + 0{A) . 



(3.9) 



[Indeed, we have (p^ = 0(3).] Similarly, by inserting (p.9| ) and ( p.31|) into the Oi component 
of (B.5|), we obtain 



Vr' = V,-jdi^' + 0{2) 



(3.10) 



Note that (|3.9| ) is valid with the inclusion of relativistic corrections e"^, while ( p.lO|) is valid 



at the non-relativistic level only. 

The relations (|3.9|) and (|3.10|) are numerically true in the region Di (1 D^. We now 
transform them into matching equations, i.e. equations relating mathematical expressions 
of the same nature. To do this, we need only to replace the inner potentials in the right- 
hand-sides of ( p.9| ) and ( |3.10| ) by their multipole expansions valid outside the source. This 



is simple because V and Vi are the retarded integrals of the compact-supported mass and 
current densities a and cTj (see ([^.3D). Thus the multipole expansions J^iV) and J^iVi) of 
V and Vi are given by 



oo 





(- 


r Y 


.r 




c) . 




(- 


r Y 


r 




c) . 



(3.11a) 
(3.11b) 



where the multipole moments V (t) and V^if) are given by explicit integrals extending over 
the mass and current densities in the source, namely 



V^(t) = / d'yvL / ^dz5,{z)a \y,t + z ^ 



V^it) = J d^yvL J^^dz6e{z)a, |^y,t + zM^ . 

In (|3.12|) , {jL denotes the tracefree part of i/l = UiiUii ■ ■ -yie, and Se{z) is given by 
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(3.12a) 
(3.12b) 



Hz) = ^^{l-z'y ; f_^dz6e{z) = l. (3.13) 



The formulas ( p.ll| )- (p.l3| ) have been proved in the appendix B of [0]. Now, we have V 



A4(y) and Vi = Aiiyi) in Di fl De, hence we can write 

= M{V) + cdtv^ + 0(4) , (3.14) 
V^-t = A<(V^,)-^5i</.o + 0(2) . (3.15) 

These matching equations relate the multipole expanded external potentials V^^^*, V^°^* given 
by ( p.l8|) to the multipole expansions M.{y), M-iVi) of the inner potentials, and can be 



used to obtain the expressions of the "canonical" moments M^lt) and Silt) entering (|2.18|) 
in terms of the source's parameters, with first relativistic accuracy in Milt) and non- 
relativistic accuracy in Silt). One can also obtain in this way the coordinate change <y9° (in 
the form of a multipole expansion). This was the method followed in [|l|, and we can check 
again that the results of |1]] are indeed equivalent to ( p.l4| )-( prT5D . (We shall in particular 
recover below the results of by a more general method.) 

C. Matching of the non-compact-supported potential Wij 

By inserting in the spatial components ij of the transformation law ( p.5|) the inner and 
outer fields (|2l9| ) and ( p.31|) , we obtain 



= W,, - j[d,^' + d,^ - 5,,(9o</^° + d,^^)] + 0(2) . (3.16) 

This equation is valid at the non-relativistic level only. Like ( p. 91) and (|3.1CI|) , it is numerically 
true in the region (1 Df.. To transform ( |3.16| ) into a matching equation, we must first 
compute the multipole expansion, valid outside the source, of the inner potential W^j which 
is, contrarily to the potentials V and Vi, of non-compact support. 



Recall that W^J^^ in the left-hand-side of ( p. 16] ) is the sum of the potential V^J"^ parametriz 



ing the linear metric ( p2.17cD , and of the nonlinear potential Z^^^ = Zij{V^^^) defined in 
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( p.29| ). Now V^*^^* has been matched to the multipole expansion of V in the previous subsec- 
tion: V^^^ = J^{V) + 0{2) as deduced from ( p.l4| ). Thus we can write the external potential 



W,f^ as 



p^ext _ ye.t ^ Zi,{M{V)) + 0(2) 



(3.17) 



where Zij{A4{V)) is given by 



Z,,{M{V)) = FPb=o Qi?' 



-d,MiV)d,MiV) + - 6,,dkMiV)dkMiV) 



(3.18) 



Note that it is crucial to replace in ZijiV^^^) not by V but by the multipole expansion 
M.(y) of V (modulo 0(2)). Indeed, contrarily to = \^ + 0(2) which is true only in the 
region Di fl Dg, the matching equation V^^^ = A4(y) + 0(2) is an identity which is valid 
"everywhere" and in particular on the whole past null cone, issued from the considered field 
point, on which depends the retarded integral in ( ^.18] ). 

As for the inner potential Wtj, we recall that it is given by 



W.J = n-^' 



-AnGai, - diVdjV + ^5i,dkVdkV 



(3.19) 



We shall now prove that the multipole expansion M.{Wij), valid outside of source, of the 
(non-compact-supported) potential ( p.l9|) reads as 

■1. 



MiW,,) = Z,,{M{V)) + V 5l 



r \ c 



(3.20) 



where the first term is given by 
second term are given by 



and where the multipole moments Wj^(t) in the 



a., + ^ (diVd.V - ^6,,jdkVdkV 



y,t + z- 



c 

(3.21) 



where 6e{z) is defined in ( |3.13| ). 

The proof of ( ^.2U| )-( p.21| ) goes as follows. The first term in Wij, involving the compact- 
supported matter stresses aij, can be treated by formulas identical to (|3.11| )- (|3.13|) . We thus 
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focus the proof on the second term involving diVdjV — from which we easily deduce the 
last term. Let us consider the difference between (minus) this term and the corresponding 
term in ( |3.18| ) involving the finite part at i? = of the retarded integral, i.e. 

X,, = n^'[d,Vd,V] - FFb=o D],^[r''dM{V)djM{V)] . (3.22) 

The analytic continuation factor in the second retarded integral deals with the singular 
behavior of the multipole expansions near the spatial origin r = 0. The first integral does not 
need any regularization factor because the integrand in perfectly regular at r = 0. However, 
in order to make a better comparison between the two integrals, let us introduce the same 
factor r^, and the "finite part" prescription, into the first integral. Since this integral is 
convergent, the "finite part" prescription simply gives back the value of the integral. Hence 
we can rewrite ( |3.22| ) as 

X,, = FFb=o □^i[r^(9,V5,V - dM{V)d,M{V))] . (3.23) 



The important point is that, under the form (|3.23|) , we see that Xij does have a compact 



support limited to the distribution of matter in the source. Indeed, outside the compact 
support of the source (i.e., in the domain Dg), the potential V numerically agrees with its 
multipole expansion Ai{V) and therefore the integrand in (|3.23| ) is identically zero. (This 
is at this point that we need to assume that the multipole expansions in the "canonical" 
construction of the exterior metric in Dg involve an infinite number of multipoles.) We can 
thus compute the multipole expansion Ai{Xij) of ( |3.23D in D^. by exactly the same formulas 
( p.ll| )-( prT^ ) as was used for the compact-supported potentials V and Vi. The only difference 
is that the analytic continuation factor must be kept inside the integral, and that one 
must apply to the formulas the "finite part" prescription. One is thus led to the multipole 
expansion 



-xh t-- 



r '"-^ \ c 



(3.24) 



where the multipole moments X^-{t) are given by 
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(3.25) 



[The "finite part" commutes with the derivation operator Ol in ( p.24| ).] We now prove that 
the second term in ( |3.25| ), which involves the multipole expansion M.{y) in the integrand, 
is in fact zero by analytic continuation. Indeed, the multipole expansion M.(y){y ^t) is by 
( p.llal ) an expansion of the type SnL/|y|~^jF(t — |y|/c), where V is some multi- index, p is 
some integer and some function of time (all indices suppressed). Thus, M.{y){y, t+z\y\/c) 
is of the type T,nLi\y\^^J^{t + (z — l)|y|/c). By differentiating and squaring the latter 
expansion, we see that the product diJ^{y)djM.{y){y,t + z\y\/c) entering ( p.25| ) is an 



expansion of the same type SnL//|y|^''^(t+ (2; — l)|y|/c), with L" a multi- index, q an integer 
and Q a function of time. We then expand by Taylor's formula each function Q when 
c — > +00. This introduces many powers of |y| and of (^ — 1). Multiplying by 5i{z) and 
integrating over z, we find a multipole expansion of the type Srij^// |y|'^7i(t), where k is an 
integer and Ti another function of time; then multiplying by / d?y\y\^yL and performing 
the angular integration (using d?y = d\y\\y\'^dQ and / dVLfiLnfiL = const. Se^in), we arrive at 
a multipole expansion of the type Tjhi,'H{t) /o^°° c^lyl lyl^"*"™, where m is an integer. Finally 
each of the latter integrals /q^°° c^lyl |y|'^'*'™ is zero by analytic continuation. [Indeed, we cut 
the integrals into two pieces, Ji = /q^ c^lyl lyl^"^™ and I2 = /y °° c^|y| |y|'^^'") where Y is some 
constant > 0. By choosing the real part of B to be such that ReB + m > — 1, we compute 
Ji = y-B+m+iy' (^^_|_^_|_]^-). Qj^^ choosing the real part of B to be such that ReB+m < —1, 
we compute I2 = — y-^+^+i/ (B + m + 1). Then both Ji and I2 can be analytically continued 
for all complex values of B except the single value — m — 1. The integral /o^°° (ilyl ly]^^™' 
is the sum of the analytic continuations of Ji and I2, and is thus identically zero on the 
whole complex plane (including —m — 1).] Note that the proof that the second term in 
( p.25| ) is zero can be easily extended to the case where we have, instead of diJ^{V)djA4(y), 
an arbitrary nonlinear multipolar product of the type d^J^idi,M.2M.3 ■ ■ ■, where J^i{y,t), 
■M.2{y, t), Ai^iy, t) ■ ■ ■ denote formal multipolar and post-Newtonian expansions of the type 
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SnL'|y|^(ln |y|)'^/C(i), where p, q are integers (such expansions are known to arise in higher- 
order nonhnear approximations of the external field [^). Thus, the multipole moments 
A'j^(t) in ( p.25| ) can be simply written as 

A;5(t) = ^ FPb=o / d'ylyfvL j\z5,{z){d,Vd,V) fy, t + zM^j , (3.26) 



and the expressions ( p.20D -( p^2l|) are then easily deduced. Indeed, we write Wij as the sum of 



Zij{A4{V)), of the compact-support integral 0~^^[—ATTGaij] and of the term —Xij + ^SijXkk- 
The multipole expansion A4(Wij) is thus the sum of Zij{Ai(y)), which is already a multipole 
expansion, of the multipole expansion of [— 47rG'crjj] computed by formulas like (|3.11|) - 
( CT ), and of -M{X,j) + ^6ijM{Xkk) as deduced from and (ICTI) . We add the factor 



|y|^ and the "finite part" prescription into the moments of the expansion of Oj^^l—AnGa, 



(which are convergent anyway), and the result is J^(Wij) given by (|3.20|) - (|3.21|) . 

One should note the remarkable role of the analytic continuation factor lyl"^ in (|3.25|) and 
( p.26| ). The integral ( p. 251 ) is of compact support and is thus perfectly well-defined at infinity 



y| +00. The role of the factor |y|-^ in (|3.25D is to deal with the singular behavior at the 



origin |y| = of the second term in the integrand, involving the multipole expansion A4{V). 
That is, one can compute ( |3.25| ) by choosing Re _B to be a large positive number so as to kill 
the bad behavior of the integral at the origin, then analytically continuing the integral near 



B = and deducing the finite part at i? = 0. On the contrary, the integral (|3.26| ) is perfectly 
well-defined near the origin |y| = 0, but it is not apparently of compact support. The factor 
lyl"^ in ( |3.26|) then deals with the a priori bad behavior of the integral at infinity |y| — > 00, 



where jjLdiVdjV behaves for large i like a large power ~ |y|^ ^ blowing up at infinity. That 
is, one can compute (|3.26| ) by choosing Re -B to be a large negative number so as to make the 



integral convergent when |y| —>■ 00. These two different procedures, ReB a large positive 
number in ( |3.25| ) and a large negative number in (|3.26| ), give the same numerical result, as 



we have just proved. In conclusion, the form ( |3.26| ) although not apparently of compact 
support is however, thanks to the properties of analytic continuation, numerically equal to 
the compact-supported form ( p.25|) . It is evident that the form (|3.26| ), where the analytic 
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continuation factor deals with the behavior of the integral at infinity from the source, is the 
one which should be used in applications. 

The equation ( |3.16| ) can finally be replaced by the matching equation 

^4 



M{W,,) - j[d,if' + d,^' - 5,,(9o¥^° + dk^")] + 0(2) 



(3.27) 



where the multipole expansion A4{Wij) is given by (|3.20| )- (|3.21| ). This equation can be 
equivalently rewritten as 



£=0 



dr 



-W^ (t - 



4 



[d,^' + dj^' - (5i,-(9o¥^° + dkv'')] + 0(2) , (3.28) 



and in this form can be used to compute, if desired, the vector (since ip^ is known from 
( p.lSp ). To this end, one needs to decompose the moments W^^ into irreducible STF tensors 
with respect to the indices ij and L. 



D. Matching equation at the post-Newtonian order 



We now have in hand all the material needed to match the inner metric ( |2.14| ), namely 



167rG 



+ 0(8,7,8) 



(3.29) 



where A and A"^ are defined in ( |2.11 ) and ( |2.12| ), to the corresponding outer metric ( ^.34 ), 
namely 

= GhZmi^cn) + FPs=o \r^T\V^^\ Vr\ W-^')] + 0(7, 7, 7) , (3.30) 



where /ican{i) linear metric (|2.17|) - (|2.18|) and A is the same expression as in ( |3.29|) but 

expressed in terms of the external potentials V^''^*^ Vf^^ and W^j^^. The metrics ( p.29|) and 
( p.30| ) are given in their respective coordinate systems and x'^^^, and are valid respectively 
in Di and DiCi D^,. 

In the two previous subsections, we have related the external potentials V^^^, Vf^^ and 
W^^^ to the multipole expansions A4(y), Aiiyi) and J^(Wij) of the internal potentials V, 
Vi and Wij. These relations, which are ( p.l4|) -( p35|) and ( p.27|) , allow us to compute the 
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effective nonlinear source A'^'^{V^^^,Vf^^,W^^^) in the right-hand-side of the outer metric 



( p.30| ) in terms of the inner potentials. Indeed, by using ( p.l4|) -( pyT5l) and ( p. 27] ) into the 
expression (gl2D of A''^ and by using DM{V) = aM{Vi) = and D*^^ = 0(7,8) (see 
( CT) ), we find 



(3.31) 



where in the right-hand-side extra terms appear which are due to the coordinate transfor- 
mation between the inner and outer metrics, and which can be written as a d'Alembertian 
operator □ acting on the tensor 



c 

Q^' = , 
Q'^ = . 



8 



(3.32a) 
(3.32b) 
(3.32c) 



[The only non-zero component is which is of order 0(6).] Now, to evaluate (|3.30| ) we 
must apply on both sides of (|3.31| ) the regularized inverse d'Alembertian operator FPn]j^. 
One easily checks that FPn]^^(nfi°'') = + 0(8). Indeed, this follows from the structure 
of which is made at leading level of terms of the type S?7,ir^(^+^+^'^\ where is a 
positive integer (see |T|] p. 395). [M.{V), M.{Vi) and ^9°, yj* all have the structure En/,/r~*^^'+^^ 



at leading level.] Hence, we find that the outer metric (|3.30|) reads as 



^canl-^can) 



Gh 



IJ,U 

caii(l) 



+ FP 



B=0 



r''A^-{M{V),M{V^),M{W,,))\ + + 0(7, 7, 7) 

(3.33) 



On the other hand, we have shown in 5 III A that the inner and outer metrics should be linked 



by a coordinate transformation involving nonlinear terms at order in its 00 component, 
see (pD. Using = -^M{V) + -- ■ and h^' = -^M{Vi) + -- ■ as is appropriate in D.nDe, 
it is easy to see that the nonlinear terms in (|3.8| ) are precisely equal to given by ( |3.32aD , 
modulo 0(8). Hence, the coordinate transformation ( p.5|) -( pl8D is in fact given by 



h^-^{x) = /i^'^(x) + d^'^" + n^'' + 0(8, 7, 8) 



(3.34) 
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(where the hnear part dip^" of the coordinate transformation is defined by ( |3.6|) , and where 
both sides are expressed in terms of the inner coordinates x^). Substituting in the left-hand- 
side of ( |3.34| ) the outer metric ( |3.33| ), and in the right-hand-side the inner metric ( p.29| ), we 
arrive at the following equation for the external linear metric (in the coordinates x^) : 



'\P,nYf_ 



-FPb=o n-^\r^A_^\M{V),Mm,M{W,,))] + d^^'' + 0{7JJ) . (3.35) 

The nonlinear part Vt^'^ of the coordinate transformation has been cancelled, and it remains 
only the linear coordinate transformation dip'^'^. 

The equation ( |3.35| ) is numerically valid in Di fl D^, but must now be transformed into 



a matching equation. The reasoning is exactly the same as the one we followed in dealing 
with the non-compact-supported potential Wij in § [III C| , and we simply repeat here the 
arguments. First, the term in ( |3.35| ) involving the compact-supported matter stress-energy 
tensor T'^'^ is treated by formulas such as ( |3.11| )-( p.l3| ). Then consider, analogously to ( p. 221 ), 



the difference between the two terms involving A^'^ in ( |3.35| ). By adding the factor and the 
"finite part" prescription in the first of these terms (whose integrand is regular at r = 0) we 
can write this difference, analogously to ( |3.23D , in a manifestly compact-supported form. (In- 
deed, the potentials V, Vi and Wij are numerically equal in De to their multipole expansions 



M.{V), M.{Vi) and A4{Wij).) Hence, we can also apply to this difference the formulas (|3.11 



( p.l3| ) and obtain, analogously to (|3.24|) - (|3.25|) , some expressions for the multipole moments 



allowing for the analytic continuation factor lyl"^ and the finite part at i? = 0. Finally the 
contribution associated with the multipole expanded source A'^'^ {Ai{V), M.{Vi), M.{Wij)) in 
these multipole moments is shown, analogously to ( p.26| ), to be zero by analytic continua- 
tion. This follows from the fact that A'^^MiV), MiVi), MiWij)) is a sum of quadratic or 
cubic terms dpAiida-J^2 or dpAiida-J^2-M.3 where the A^^'s admit formal multipolar and 
post-Newtonian expansions of the type Sr;,i^/|y|^(ln |y|)^/C(t), with p, q some integers {q = 
in the case of A4(y) and J^iVi), but g = or 1 in the case of A^(VFjj)); see the discussion 
above (|3.26| ). An important result of this paper can now be written down. We introduce an 
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effective (truncated) total stress- energy tensor of the matter fields and of the gravitational 
field, 



iOTrCr 

where A and A.^" are given by ( p.ll| ) and ( |2.12| ). By ( |2.13| ) this tensor satisfies 

9,r'^^ = 0(3,4) . 



(3.36) 



(3.37) 



Then we can transform ( p.35| ) into the matching equation 



G'C,i,|A^L.Stl 



^=0 



-TTit-'- 



+ dip^' + 0{1, 7, 7) 



(3.38) 



relating the exterior linear metric (|2.171 )- p.l8| ), that we recall is a functional of the "canoni- 
cal" moments Ml and S'l, to the total truncated stress-energy tensor ( p.36| ) via the multipole 
moments 



(3.39) 



The result ( p3.38D -( P75UD is especially simple. It says that the linear metric /ican(i) 
equal, modulo the linear coordinate transformation d^p^'^ (and modulo 0(7, 7, 7) terms), to 
the multipole expansion outside the source we would obtain z/the effective total stress-energy 
tensor T^^ had a compact support limited to the material source only. The difference is that 
the moments ( |3.39| ) carry the analytic continuation factor |y|'^ and the "finite part" at 5 = 
to deal with the poles at 5 = coming from the behavior of the integral at its upper bound 
|y| +00. When no poles arise as will be the case at the 2-PN approximation (see § [TVl ), 
we can say that p.38| )-( |3l3^ ) justifies the formal procedure followed by Epstein and Wagoner 
30| and Thorne to compute the multipole moments. (However recall that in the works 



|30| and |^ these multipole moments are assumed to be the moments which are radiated at 



infinity, while one has still to add to these moments all tails and nonlinear contributions in 
the radiation field, see § [TV] below.) Let us emphasize again that the expression ( |3.39| ) is not 
manifestly of compact-supported form but is numerically equal to the expression obtained 
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by multipole expanding the compact-supported right-hand-side of ( p.35| ). In particular, this 
shows that the multipole moments ( |3.39| ) are retarded fnnctionals of the source's parameters, 



depending on the source at times t' < t only (or t' < t — r/c in ( |3.38| )), contrarily to what 
is apparent on their expressions ( p.39| ) but in accordance with what they must be. 



IV. GENERATION OF GRAVITATIONAL WAVES 



A. Relations between the canonical moments and the source moments 



In order to find from the matching equation (|3.38|) -( ^.39 ) the expressions of the "canoni- 



cal" moments M^, Sl as functions of the source's parameters, it suffices simply to decompose 
the reducible moments T'^(t) into irreducible STF multipole moments. This has already 
been done by Damour and Iyer [^| in the case of linearized gravity, i.e. in the case where the 



total stress-energy tensor t^'-' in (|3.38|) is replaced by the usual compact-supported stress- 



energy tensor T^'^ of the matter fields, supposed to be exactly conserved: duT^^" = 0, and 
where of course there are no analytic continuation factors in the expressions of the moments. 



Let us prove that the computation done in |31| basically applies to our case, which 
involves both analytic continuation factors and a stress-energy tensor r^^" which is only 
approximately conserved: d^^T^^'^ = 0(3, 4) by ( ^.37] ). To this end, we need only to check that 
the multipole expansion in the right-hand-side of ( |3.38| ) is divergenceless up to 0(3,4) as 
a consequence of the approximate conservation of T'^". We take the time derivative of the 



0/i components of the multipole moments ( p. 391) , use c ^dtT^^ = —djT^^ + 0(3,4), perform 



some integrations by parts both with respect to y and to z, and finally insert the identity 

{d/dzy6e-^-i{z) = (2£ -|- 1)(2£ -|- 3)[6e-i{z) — 6e{z)]. These operations result in 



2 



+FPb=o[b J d'y\yf-%yL f^dz6,{z)T^^^{y,t + z\y\/c)^ 
+0(3,4). (4.1) 
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The first two terms in the right-hand-side of ( [4.1| ) are the ones which would ensure the 
exact zero divergency of the multipole expansion in ( |3.38| ). {Tj^^l^^ means the STF part of 
Tl^l\-) The third term is the finite part at i? = of an integral having explicitly i? as a 
factor. This term arises because of the derivation of the factor lyl"^ during the integration 
by parts. (Note that we have discarded a surface term which is easily seen to be zero by 
analytic continuation.) Finally, the remainder 0(3,4) comes from the only approximate 
conservation of T'^'^. We now show that the third term in ( [4. 1]) in fact vanishes up to order 
0(3,4), which means, thanks to the explicit factor B, that the integral has at this order no 
simple pole a.t B = 0. This follows from a general lemma which we now state and which 
will turn out to be very useful in the remaining of the paper. Let y be a "source" point, xi, 
X2,. . . , Xn be n "field" points, and ao, ai,. . . ,an be n -|- 1 real numbers. Then we can state 
that 

FPb=o [b J d'y\y\''+'''yL\y - xil"^ ■ ■ • |y - x,r"| = (4.2) 

if the sum ao + ai + ■ ■ ■ + is not an odd (positive or negative) integer. To prove that 
this is true, we need to investigate the behavior of the integrand at the bound |y| — >■ 
+00 of the integral, where each |y — xil"^,. . . admits a multipole expansion of the type 
^2;^i-i^^jy|ai ... (coefficients suppressed). Thus, one is led to consider the behavior of 
integrals of the type / d^y\y\^~^'^°yLdLi |y|"^ ■ ■ ■ •^Lnlyl"'"- By performing the integration over 
the angles (using d^y = \y\'^d\y\dfl), we find that these integrals are zero if i + £i + ■ — h 
is an odd integer. When i + ii + ■ ■ ■ + £^ is an even integer, a single type of radial integral 
remains, which is / d\y\ \y\^'^^ where P = 2 + ao + £+ (ai — ii) + ■ ■ ■ + (on — ^n)- A pole 
at i? = will arise only if (3 = —1, and thus, since £ + ii + ■ ■ ■ + is an even integer, only 
if ao + en + ■ ■ ■ + ctn is an odd integer. If ao + ai + ■ ■ ■ + an is not an odd integer, there 
is no pole and ( [4.2| ) is true by virtue of the explicit factor B in front of the integral. (Note 
that this condition for ([4. 21) to hold is sufficient but by no means necessary.) Now, a non- 



zero contribution to the third term in (|4.1|) can come only from the non-compact supported 



part of T^^ in the integrand, i.e. the part involving A""^ (see ( p.36|) ). We first expand the 
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argument t + z\y\/c when c +00 up to an order consistent with the remainder 0(3, 4) in 
( [4.1[ ), and integrate over z (using ( [4 .81 ) below). Thus we have to consider integrals involving 
A''^ and its time- derivatives, and multiplied by some even powers of |y|. The structure of 
A''^ as given by ( |2.12b| ) and ( |2.12c| ) is of the type dVdV, where V represents a retarded 



potential ( |2.3a| ) or ( |2.3b| ) and d is some time or space derivation (all indices suppressed). By 
expanding the retardation argument in the V^'s up to the same order 0(3, 4), we see that the 
structure of A""^ is of the type dUdU or dUdX, where U and X are instantaneous (Poisson- 
type) potentials of the type / (i^xcr(x)|y — x|" with a = ±1 (see ( [4.14| )- (|4.15| ) below). Thus, 
A"*^ is composed of terms of the type / cP:x.id^X2dx-i^dx20'{xi)(7{x2)\y — xi|"i|y — x.2\°'^ where 
ai + 0^2 = —2 or 0, and where the derivatives act on the source points xi = (xi,t) and 
X2 = (x2,t) (we use the fact that (9y|y — x|" = — 9x|y — x|"). Finally, it remains to commute 
the integration over y with the integrations over Xi and X2 to arrive at a series of integrals 
of the type / d^'y\'y\^~^°"^yL' |y |y — X2|"^, where ao is an even integer and ai + 02 = —2 

or (we have used i/iyL ~ yn + |yP^j<jf yL-i>)- These integrals have no poles by the lemma 
]2|). Therefore, we have proved that (|4.1|) reads in fact as 



±nit) ^ .TS, W + ^ A T-«) + 0(3. 4) . (4.3) 



This relation shows that the multipole expansion in the right-hand-side of ( p.38|) is diver 



genceless modulo small terms of order 0(7,8). This is consistent with the fact that the 



"g-part" of the external metric is found to be zero at order 0(7,7,7), see (|2.25|) and the 
appendix 0. Because of the uncontrolled remainder in ([4.3|) , we find that the four relations 
(5.27) in are modified by small post-Newtonian terms: Cl is modified by a term 0(7), 
and Ql, Tl, Jl are modified by terms 0(8). Thus the Oi components (5.9b) in involve 
some uncontrolled terms 0(7), while the ij components involve uncontrolled terms 0(8) (the 
00 component not being affected). All these terms fall into the remainder term 0(7, 7, 7) in 
( p.38| ), showing that the end formulas (5.33) and (5.35) of giving the expressions of the 
mass-type and current-type multipole moments in linearized gravity can be applied to our 
case with the replacement of the matter stress-energy tensor T^^ by the total stress-energy 
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tensor T'^", and with the addition of analytic continuation factors |y|^ and of the "finite 
part" prescription. 

The matching equation ( |3.38| ) can thus be written in the form 

GKL(i) [Ml, S,] = GKL,, [h, Jl] + du^'^ + d^^'^ + 0(7, 7, 7) , (4.4) 



where duo^^ is a certain harmonic hnear gauge transformation (□cu'^ = 0), and where the 
mass-type and current-type source moments II and Jl are given by 

4(2£+ 1) 



hit) = fPb=o I t^Vly"" 



dz 



2(2£+ 1) 



c2(£+l)(2£ + 3) 



y,t + z- 



(4.5a) 



j^it) = fPb=o / d'y\y 



iB 



dz 



21+1 



c2(£ + 2)(2£ + 3) 



^l{z)£ab<iiyL-l>a^b 

^i+l{z)£ab<iiyL~l>acdt^ 



be 



y,t + z- 



The mass, current and stress densities in (|4.5| ) are defined by 



y.. —T^J 



(4.5b) 



(4.6a) 

(4.6b) 
(4.6c) 



where the stress-energy tensor r^"^ is given by ( |3.36|) . The weighting function 5i{z) is defined 
in ( |3.13| ). The matching equation ( [4.4| ) then tells us that the gauge transformation duj^" 
necessarily satisfies duj^^ = —dip^^ + 0(7,8) (indeed remember that □y^'^ = 0(7,8)), and 
that the "canonical" moments Ml, Sl in the left-hand-side of ( [4.4| ) are related to the source 
moments II, Jl by 



Mi(t)=/i(t)+0(5) , 

SL{t) = JL{t)+0{4). 



(4.7a) 
(4.7b) 
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[The relation ( |4.7a| ) for the mass moment comes from the 00 component of the matching 
equation ( |4.4| ), while the relation ( [4.7b| ) for the current moment comes from the Oi compo- 
nents of ( [4.4[ ).] The relations ( [4.7a| ) and ( [4.7b| ) are exactly the ones needed to solve the wave 
generation problem at the second-post-Newtonian approximation. 



B. Expressions of the mass-type and current-type source moments 



Since the relations ( [4 .71) are only valid up to some post-Newtonian order, it is sufficient 
to consider the post-Newtonian expansion c +oo of the source moments I^, Jl defined 
by (|4.5| ) . This expansion is achieved by means of a formula derived in the appendix B of [|l| , 
and giving the expansion when c +oo of terms involving the average over z appearing in 
rsl). This formula reads as 



dz5e(z)Il y, t + z — 
1 V c 



S(y,t) + 



2c2(2£ + 3) 



8c4(2£ + 3)(2£ + 5) 
Using (|4.8| ), and retaining consistently the powers of c^^, we obtain 



9,^S(y,t) + 0(6) 



(4.8) 



hit) = FFb=o I d'y\y 
4 



B 



\y?yL 



2c2(2£ + 3) * 



lyl^^L 



k4(2£ + 3)(2£ + 5) 



2(2£+l)|y|2y, 



iL 



c4(£ + l)(2£ + 3)(2£ + 5) 



c2(£ + l)(2£ + 3) 



c4(£+l)(£ + 2)(2£ + 5) 



(y,t) + 0(6), (4.9) 



JL{t) = FPB=o£ab<ii / d^y\y 



iB 



v=r I y yL~l>a o2vr 



't^bc 



2c2(2£ + 3) 

{y,t) + o{A) 



(4.10) 



c2(£ + 2)(2£ + 3) 

We must then insert into ( [4.9D -( ^!T0l) explicit formulas of S, Ej and Sj^ which are easily 
computed from ( p.36[ ) where we use the expression ( p.ll| ) of A = |(7|+0(6) and the expressions 
( p.l2| ) of the components of the effective nonlinear source A^" . We find the formulas 
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1 + 



diVd^y 



^ / - V9,V - 2VidtdiV - WiMV - i (dtvy 



GcH 



7. 



+29iVj9jri + 2diVdiW - -VdiVdiV 



1 + 



AV 



a. + ^ [duV{d,Vk - duVi) + \ + 0(4) 



S., = a,, + ^ ja.yo'.l^ - ]^6ijdkVdkV^ + 0(2) . 

Again retaining the powers of consistently with the accuracy indicated in 
we can now write formulas for the source moments in raw form: 



(4.11a) 
(4.11b) 
(4.11c) 

D-O, 



hit) = fPb=o / rfV|y| 



B 



i + -v--{w-v' 



1 



ttGc^ 



VLdiVd^V 



TT 



Gc^ 



Vl 



-Vd^V - 2V,mV - Wi.df^V - -{dtVy 

+2diVjdjVi + 2diVdiW - IvdiVdiV 

2 



\y?yL 



2c2(2£ + 3) 



1 + — \ a 



TdiVd,V 



\y\^yL 



8c4(2£ + 3)(2£ + 5) 
Ai2i+mL 



dt 



c2(£ + l)(2£ + 3) 

2{2i+l)yi^L 
'c^{e + l){£ + 2){2i + 5) 



2i2i+l)\y\^y,L 
c4(£+ l)(2£ + 3)(2£ + 5) 
1 

^. \ rr- -I 



t ^« 



. 4V , 

l + — ](T. 

c 



^ 2 1 dkV{diVk - duVi) + -dtVd.V 
ttGc^ I 4 



0(6) 



(4.12) 



|y| yL-l>a q2 



2c2(2£ + 3) 



vr 



Gc2 



yL-l>a 



dkV{d,Vk - dkVh) + -dtVd,V 



{2i+l)yL-l>ac ^ 

"c2(£ + 2)(2£ + 3) * 



47rG 



dbVd.V 



0(4). 



(4.13) 



Recall that the analytic continuation procedure in ( |4.12| )-( P:.13| ) is needed only for handling 
the (apparently) non-compact supported terms like the term yLdiVdiV in ( [4.12| ), and could 
be removed from manifestly compact supported terms like e.g. yicrV'^. 

Let us now derive an equivalent but somewhat simpler form for the mass-type multipole 
source moment (H:-12|) . First we replace the retarded potentials V, Vi, Wij used in ( [4.12|) 
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by "instantaneous" potentials U, X, Ui, Pij (and P = Pa) defined by the post-Newtonian 
expansions 

1 



V = U + —d^X + 0{3) 
Vi = U, + 0(2) , 
W,i=Pi,+0{l) . 



These instantaneous potentials are 

^(x,t)=G//^^a(x',t) 
J X — x' 



X(x, t)=G J t/V|x - xV(x', t) 

|x — x'l 
rf^x' 



P,,{^,t)=Gj^^_^ 



(^ij + {diUdjU - ^SijdkUdkU 



(x',t) , 



X — X 



era - ^crU 



. . ^ 



We can then rewrite ( [4.12| ) as 
hit) = FPb=o J d'y\y\ 



1 + ^f/ + ^{2d^X -8P + 8f/2) 



-diUdid^X - Ud'fU 



-2Uidtd,U - P^MU - -{dtUf + 2diUjdjU^ + 2dMP - -UdiUdiU 



\y\ vl ^2 

2c2(2£ + 3) * 



1 + 



4U 



; ttGc^ 



diUdiU 



_5V 2(2£ + l)|yP^., 

8c4(2£ + 3)(2£ + 5) * c^{i + l){2i + 3){2i + 5) * ^ 



'iL 



c2(£+l)(2£ + 3) 

2{2^ + l)m,L 



1 



ttGc^ 



l + ^]ai + { dkUid^Uk - dkU{) + ^dtUdiU] 



(4.14a) 
(4.14b) 
(4.14c) 



(4.15a) 
(4.15b) 
(4.15c) 

(4.15d) 

(4.15e) 



+ 0(5). 



(4.16) 



c4(£+l)(£ + 2)(2£ + 5) 

Several reductions of this expression can be done by integrating various terms by parts. 
For instance we can operate by part the term ~ yidiU diU using /\yL = and AU = —AnGa. 
This yields 
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FPb=o / d^y\y\''yLdiUdiU = YVb=, J d'y\y\'' [^nGyLaU + ^9, [yLd^U^ - d.yLU^] | . 

(4.17) 

The second term in the right-hand-side of ( [4.171 ), which is made of the product of |y|^ and 
of a pure divergence, is easily shown to be zero thanks to the lemma ( [4 .21) . Indeed, the 
differentiation of the factor |y|^ yields = —2, and since the term involves two potentials 
U it has «! + ^2 = —2. Hence we can write 

FPb=o J d'ylyl'^yLd^UdiU = 4nG J d'yy^aU , (4.18) 

where we have removed the factor |y|^ in the right-hand-side since the term is compact- 



supported. Note that the use of the lemma (|4.2| ) permits to freely integrate by parts all 



non-compact-supported terms in ( [4.16| ) except the term involving the product of three t/'s, 
i.e. the term ~ yiJJdiUdiU , which must be treated separately. For this term we write, like 
in (|A3), 

FPb=o / d'y\y\''yLUdiUdiU = FPb=o / rfVlyl"" [^nGyLaU^ + ^8, [yLd^U'' - d.yLU^] | . 

(4.19) 

To cancel the second term in the right-hand-side, we have to show that the integral 
/ '^^y|y|'^~^yL'|y ~ ^il^^ |y ~ ^2!"^ |y — ^3!^^ has no pole at 5 = (this case is not cov- 
ered by ([4. 2D). When |y| +00 the integrand behaves like \y\^^^'~^^~^'^~^^~^hLihL^nL^hL^ 
(where rii = yi/\y\). The angular integration shows that ii + £2 + ^3 = ^' + 2p where p is a 
positive integer, so that the remaining radial integral is / c?|y| lyl"^"^^"^ which cannot have a 
pole (this would mean p = —1). Thus we can conclude, like in (f4.18|) , that 



FPb=o / d^ylyl^'yLUdiUd.U = 2ttG J d^y^aU^ . (4.20) 

Incidentally, the above proof shows that in higher post-Newtonian approximations, terms 
can arise which cannot be integrated by parts without generating poles. This is for instance 
the case of the term ~ yiJJdiUdid^X which is expected to arise at the 0(6) level. 
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The identity ( [4.18| ) shows that the term ~ yidiUdiU in ( |4.16|) exactly cancels a previous 
term ~ yL<^U, so that we recover, at the first post-Newtonian approximation, the expression 
obtained in i.e. (2.27) in [Q. We use also the identity ( [4.20| ), and then perform several 
other manipulations like for instance one showing that the two terms involving the poten- 
tial X in ( |4.16| ) can be advantageously replaced by a single term ~ yiUd^U — all these 
manipulations are justified by the lemma (f4.2|) . They yield the expression 



lL{t) = FFB=o J d^y\yr<,yL 

4(2£+l)y,i 



+ ^{cTiiU - aP) 



2c2(2£ + 3) 



c2(£+l)(2£ + 3) 



1 + — ai + 



5V 2(2f + l)|yPy., 

8c4(2£ + 3){2i + 5) * c4(£ + l){2i + 3){2i + 5) * ' 



2(2£+ 1) 



(Jin H 

. AttG 



diUdjU 



c^{e + l){i + 2){2i + 5) 
^ ' ^ P^.d^U -2UAd^U + 2diUjd,Ui-hdtUf -Ud^^U 



Vl 



} + 0(5). (4.21) 



The expression ( |4.21| ) will be applied, in a forthcoming paper [Q, to the problem of the 
generation of waves by a coalescing compact binary at the 2-PN approximation. Note that 
many other transformations of the expression ( [4.2 1| ) could be done using the equations of 
motion and of conservation of mass. In the appendix |B| we show that in the case i = the 
expression ( [4.21| ) reduces to the expression of the conserved total mass M at 2-PN which is 
known from the equation of conservation of mass at 2-PN. 



Finally we prove that the expression ( [1.13| ) of the current-type source moment is equiv- 
alent to the expression obtained in We make the comparison with the expression (5.18) 
in 1^, which involves a STF tensor Y^{yi,y2) depending on two source points yi, y2 and 
defined by 

r-l 



>^nyi,y2) 



|yi -ysl / day}^ 



(L) 



(4.22a) 



where j/^ = ay{ + (1 — a)yl, and where y^^^ denotes the STF part of y^ = Ua ■ ■ -Ua- An 
alternative form of can be obtained by explicitly performing the integration over a. It 
reads 
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Y^y.,y.) = ^^Y.ytV , (4.22b) 



p=0 



where we sum over the number p of indices present on yl^ = y2''' Vi (i^i which case 
yf"^ = y^"^^ ■ ■ -yi) and without p-dependent coefficient in the sum, e.g., V^^ = ^\yi — y2\ 
^{yiUi + 2/1*1/2'* + 2/2*2/2^)- Then the formula which permits to relate our work with the 
formalism used in 101 is 



d^y\y 


IV 


|y-yi| 


y-y2| 



FPb=o / 1 , = -27rr^(yi, ya) . (4.23) 



The proof of this formula is as follows. We know from (3.9c) in that (|y — yi||y — y2|)^^ = 
AyA; — 2n6i2 where 612 = |yi — y2| Jo da6{y — y^) represents a Dirac distribution on the 
segment joining yi and y2 {6 is the usual three-dimensional Dirac distribution), and where 
k is some kernel given by A; = ^ ln[(|y — yi| + |y — y2|)^ — |yi — y2p]- By substitution into 
the left-hand-side of ( [4. 231) one gets two terms. The ffist one is the finite part at I? = 
of the integral / d^y\y\^yL^yk, which is also equal to B J d^y\y\^~^yi[diyLk — yidik]. We 
replace into the latter integral the convergent Taylor expansion when yi, y2 ^ of the 
kernel k (which is analytic in yi, y2; see (3.14) in 0), and find that the only remaining 
radial integrals are of the type B J d^y\y\^~'^'^ ot B J d^y\y\^^^ \n |y|, where k is an integer. 
These radial integrals are zero at B=0 (no poles). Thus it remains only the second term 
which is the finite part at 5 = of / d^y\y\^yL{—2n6i2) , and readily yields the result ( [4.231) . 



By multiplying ( |4.23| ) by some densities cr(yi,t) and cr{y2,t) and integrating over d^yi and 



'^^y2 we obtain such relations as 

FPb=o / d^ylyl^'yLd.Ud.U = -27^1/ d^yid'y2a{yi)a{y2)Y^^{yi,y2) , (4.24) 

where Yfj = d'^Y^ /dyldyl, showing the complete equivalence between our expression ([4.1 3|) 
above (where the ^'s can be replaced by the corresponding f/'s) and the expression (5.18) 
in m. 
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C. The asymptotic waveform at the 2-PN approximation 



It has been shown in (see also |3^) that the "canonical" external field ( |2.16|) , which 



satisfies all-over the harmonic-coordinates Einstein's equations (1J.)-(L2) (in the sense 
of formal nonlinear expansions), can be rewritten in a so-called radiative coordinate system 
= (cT, X) in which it is of the Bondi type at large distances from the source. It 
is sufficient to consider the transverse-traceless (TT) projection of the leading-order 1/R 
part of the spatial metric (where R = |X| is the distance to the source). Denoting by 



km ) 



\TT 



this TT projection of the spatial metric (where Qkm is the 



usual covariant metric), we can then uniquely decompose the 1/R part of /i^^ into the infinite 



multipole moment series 



AC °° 1 r 

hZiy., T) = ^V.,U^) E ^\Nl-2U.,l-2{T - R/c) 



2i 



NaL~2£ab{iVj)bL-2{T - R/c) 



(4.25) 



(£+l)c 

where the radiative moments and Vl represent two infinite sets of functions of the retarded 
time T — -R/c, which are STF in their indices L = ii ■ ■ - ii {£ goes from 2 up to infinity). 
These functions by definition parametrize the asymptotic waveform. The coefficients in 
( [4.25| ) have been chosen so that the moments Ul and Vl reduce, in the non-relativistic limit 
c — s> +00 to the ith time-derivatives of the usual "Newtonian" mass-type and current-type 
moments of the source pO[. Our notation in ( [4.25| ) is Nl-2 = -^ii ■ ■ ■ -^i£_2 with Ni = X^/R, 
NaL-2 = NaNL-2, T^ij) = \{Tij + Tjj), and 



V^.kmm = {5,k - N,Nk){5,m - NjNj - -{6,, - N,N^){6km - N^NJ . (4.26) 

At the 2-PN approximation, including all terms up to the level ~ c~^, the waveform 
( [4.25| ) reads as 



"'km 



2G r 1 

c^R^'"'"'[ ^ c 
1 



1 4 

-NaUija + -£ab{iVj)aNh 
~^X(if)Uij(ii) -\- ~£ ab{iV'j'jacXl)c 
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1 

+ - 

1 



QQ-^abcUijabc ~^ -^^^ ab{i^j)acd-^bcd 



r^QQ-l^abcdUijabcd ~l~ r^^^ ab{i^j)acde-^bcde 



+ 0(5) 



(4.27) 



By differentiating, squaring and averaging over angles this expression, one obtains the energy 
loss formula at the 2-PN approximation, giving the rate of decrease of the Bondi energy Eb'- 



dE 



B 



dT 



189 ■''^ •''^ 45 " 



16, 



^(i)t/(i) 



9072 "''^'" -''^"^ 34 



^(i)t/(i) 



+ 0(6) 



(4.28) 



Now the point is that the external field ( |2.16| ) is algorithmically constructed in ||21| 
from the linearized metric ( |2.17| )- p.l8| ), which is parametrized by the "canonical" multipole 
moments M^, Sl, and that the coordinate transformation between the harmonic coordinates 
and the radiative ones can also be algorithmically implemented [^,^]. Therefore, the 
radiative moments Ul and Vl parametrizing the multipole expansion of the asymptotic 
waveform ( |4.25| ) are necessarily given as some algorithmically computable functionals of the 
canonical moments Ml and Sl. It has been shown in previous papers [|l],H that Ul and Vl 
are given by some nonlinear infinite expansions in G (consistently with our whole approach) 
of the type 



n>2 



^3(n-l)+S£--£ 



XnL{T) , 



G 



n-l 



n>2 



g3(n-l)+S^.-£ 



(4.29a) 
(4.29b) 



where M^\t) and Sl^I\t) denote the ith and {i — l)th time-derivatives of Ml and Sl-i 
computed at the radiative time T, and where X„i(T) and YnL{T) are nonlinear functionals 
of order n of the moments Ml and Sl and their time-derivatives. The general structure of 
XnL and YnL is 

XnLiT),YnLiT) = J2 r dU,--- r dUn}CLL^...L^iT,U,r ■ ■ ,Un)Vi^^'\u,) ■ ■ ■Vt\Un) , 



(4.30) 
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where V^^^ denotes the ath time-derivative of either a mass moment Ml (in which case 
£ = or a current moment eai^^j^igSaL-i endowed with a Levi-Civita symbol (in which case 
£ = i + 1). The tensor ICll-^-l^ denotes some dimensionless kernel whose index structure 
is made out only of Kronecker symbols, and which depends only on variables having the 



dimension of time. The powers of G and c in ( 4.29 ) are obtained by a simple dimensional 
argument, namely that the mass and current moments Ml and 5*^ have the usual dimensions 
of multipole moments. The notation is for = + s, where is the total 

number of indices present on the n moments in ( [4.30| ), and where s is the number of 
current moments among these n moments. As the tensor ^ll^-l^ represents an operation 
of complete contraction between the indices L,L_i, - ■ ■ , L„, we have necessarily the equality 

El. = ^ + 2A;, (4.31) 

i=l 

where k is the number of contractions among the indices Li, ■ ■ -L^. 

In view of the explicit powers of in front of the multipole moment contributions 
present in ( |4.27| ), we need to compute the relations ( ^.29| ) linking the radiative and canonical 
moments only up to some definite order in c~^. Namely, Uij is to be computed up to 
inclusively, Uijk and Vij are to be computed up to c~^, Uijkm and Vijk up to c~^, and so on. 
Now the equality ([4.3 1|) shows that the nonlinear terms in (|4.29| ) having n > 3, and thus 



coming from the cubic and higher nonlinearities of Einstein's equations, are at least of order 
0(c^^) = 0(6) and can be neglected for our purpose. Furthermore, the terms with n = 2 
must have k = in (|4.31| ) since for /c > 1 the corresponding order is 0(5), also negligible for 
our purpose. Finally, the remaining nonlinear terms with n = 2 and k = 0, which represent 
corrections ~ in the radiative moments, are to be computed only in the quadrupole and 
octupole mass moments Uij and Uyk, and in the quadrupole current moment Vij. We then 
easily arrive at the only possibilities i = li + i.2 (see ( [4.31[ )) with i = 2 (case of Uij) and 
= 0, I2 = 2, or £ = 3 (cases of Uijk and Vij) and li = 0, I2 = 3 or ii = 1, £2 = 2. (Indeed, 
one of the moments in ( f4.3(J| ) is necessarily non-static, ^ 2 say.) This corresponds to the 
interaction of the mass monopole M and of the mass quadrupole Mij (case of Uij), of M 
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and of the mass octupole Mijk or of the mass dipole and of Mij (case of Uijk)-, and of M 
and the current quadrupole Sij (case of Vij). Let us combine these informations with the 
results of showing that two (and only two) types of "hereditary" contributions arise in 
the radiative moments Ul, Vl sd. the quadratic nonlinear approximation, namely the "tail" 
contributions involving the interaction between M and non-static moments or Sl, and 
the "memory" contribution involving the interaction between two non-static moments M^. 
By the previous reasoning, the latter "memory" contribution can be neglected, and the 
former "tail" contributions need to be included only in the radiative moments Uij, Uijk and 
Vij. Hence we can write, from (2.42) and (3.4) in 0, 

2GM 



U,,{T) = Mi;\T) + ^ dV 



U^AT) = Mf^l{T) + '^-^ j^^ dV 



In 






V26y 


+ K2 


In 






V26y 





M^f{T-V) + 0{5) , (4.32a) 
Mip,iT-V) + 0i5) , (4.32b) 
Sif{T-V) + 0i5). (4.32c) 



hoo 



V,(T) = Slf(T) + ?^C- dV 



C Jo 

The other radiative moments Uijkm, ■ ■ ■ , Vijkmn in ( |4.27|) and ( [4.28| ) are equal, with the re- 
quired precision, to the corresponding M^jl^, ■ • • , S^\^^. Three purely numerical constants 
K2, K3 and n'2 appear in ( |4.32| ), which are in factor of "instantaneous" (non-hereditary) contri- 
butions. Note that there is a priori also a contribution involving the interaction between the 
mass dipole Mi and the mass quadrupole Mjk(T) in the octupole moment UijkiT) of (4.32b). 
This contribution, which is necessarily "instantaneous" and of the type ((j3G'/c^)M(jMj^j(T), 
where (J3 is some numerical constant, has been set to zero in ( [4.32b| ) by requiring that the 



(harmonic) exterior coordinate system is mass-centred, i.e. Mj = 0. The computation of 
K2, and n'2 necessitates the implementation of the algorithm for the construction of the 
external metric. This was already done in for the computation of K2. The computation 
done in the appendix yields the values 

11 97 7 

Ko = — , /t'? = — , «:', = -. (4.33) 

12 60 ^6 ^ ' 



The constant h entering the tail contributions in ( |4.32| ) is a constant (with dimension of 



a time) which parametrizes the relation between the radiative coordinate system (T, R) in 
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which the metric is of the Bondi type and the harmonic coordinate system (tcan, ''"can) of 
5 lilBI. It is such that 



^ R ^ r^an '2GM fr,,^ 

T = tcan — In 1 — 



+ 0(5) 



(4.34) 



(where terms of order O^l/r"^^^) in the distance to the source are neglected). A possibly 
convenient choice for the constant 6 is 6 ~ l/tuo, where ujq is a typical frequency at which 
some detector at large distances from the source is operating ||TT|. 

The relations ( [4.32|) are still not expressed in terms of the source's parameters, and the 
last step obviously consists in using the relations (|4.7| ) linking the "canonical" moments Ml, 
Sl to the real source moments II, J^. We can thus rewrite ( [4. 321) as 



U^AT)=I\fiT) + 



2GM /•+°° 



dV 



c Jo 
2GM r+<^ 



' 26/ 12 



lif{T-V) + 0{5) 
l\%{T-V) + 0{h) 
j\f{T-V)+0{A) 



dV 



■' O" JO 

m 2GM r+°° 

V^AT) = JS\t) + ^ I dV 



In 



97 

60 
7 



2bJ 6 



(4.35a) 
(4.35b) 
(4.35c) 



(with relations limited to the first term in the right-hand-side for the higher-order moments 
Uijkm,---)- One must insert these relations, together with the explicit expressions (|4.21|) 
and ( ^4.13|) of the source moments, into the waveform ( |4.27| ) and/or the energy-loss formula 
( [4.28| ). (Note that the only tail contribution in the energy-loss formula ( [4.28[ ) comes from 



the "mass-quadrupole" tail associated with the moment Uij in ( |4.35a| ).) This solves the 
problem of the generation of gravitational waves by a general isolated system at the 2-PN 
approximation. 
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APPENDIX A: PN EXPANSION OF PART OF THE EXTERNAL FIELD 



The "g-part" of the external metric is the second term in the definition of the nonhnear 



canonical field ( 2.21 ). We shall denote the first term in ( 2.21 ) by 



(n)J 



(Al) 



The "g-part" of the metric is computed from the "p-part" (^T|) as follows [PTI] . We compute 



the divergence of (^), namely 



r 



can(n) 



d^PZn(ny and obtain 



(n)J 



(A2) 



where = x'^ /r and where the factor B comes from the derivation of the analytic contin- 
uation factor . This divergence is known to be a retarded solution of the wave equation, 
and can thus be decomposed, in a unique manner, as 



r 





can( 



£>0 



can(; 



l>0 



r \ c 
r \ c 



(A3a) 



t>i ^ 



1 / r 

-CiL-l [t 

r V c 



+ ^iabdaL- 



1 ^ / r 
-DbL-i [t-- 
r \ c 



(A3b) 



where A^, B^, Cl and are some STF tensorial functions of the retarded time. Then the 
"g-part" of the external metric is defined by its components as 



00 
ycan{n) 



ycan{n) — (^^iah<->a \ J^, 



r 



1ca,n{n) ^ij 



r 



J2dL-l(-AL- 



i>2 



-B + da (-Ba 

r \r 



4(1) , 3 ^(2) 1^ 



+ 26ijdL { ^Bl^ - 6dL-l(i - 2daL-2 (£ab{i^Dj)bL-2 



(A4a) 
(A4b) 



(A4c) 



(where e.g. A(-'\t) = t^dt'A{t')). 
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The main task is to deal with the quadratic case n = 2. We first control the post- 
Newtonian expansion of the divergence r'^^j^(^2) (0)- The quadratic source A'^!^^^2) ~ 



N^'^{hcs.n{i)) is computed by inserting the linear metric (p.l7| )-( |2?T8| ) into ( |1.5| ). Its post 



Newtonian expansion starts at 0(4, 5, 4), with a next term at 0(6, 7, 6). Let us write 

= Jt^F^' + + 0(7, 8, 7) , (A5) 

where F^'^ and are the coefficients of the leading-order and next-order terms in the 
post-Newtonian expansion, and where we use the notation u = when = 00 or ij and 
uj = 1 when fiu = Oi. Then it is easy to show that the structures of F'^'^ and G'^'^, as concerns 
their spatial dependence, are 

= dp (^) dQ (i) , (A6a) 

p+q>2 ^ ' ^ ^ ' ^ 

= Y dp (i) (i) dp (i) dc^ir) , (A6b) 

p+q>l P+<?>2 ^ ^ 

where P and Q are multi-indices with p and q indices. Important for our purpose is the fact 
that the number of spatial derivatives is p + g > 2 in F^'^ and in the second term in C^", 
and is p + g > 1 in the first term in G^"^ . Knowing the expansion (|A5|) we can write the 
corresponding expansion of the divergence ( |A2| ) . We have 



+ ^FPb=o ^A-'[Br^-'n,Gn + (^^^ A-'[Br^-'n,F'^^] }+ 0(8,7) , (A7) 

where is the Poisson operator and = (A^^)^, and where u = 1 when /z = and 
u = when fi = i. The justification of the equation ( |A7|) can be found in our previous papers 
(see e.g. (3.25) in gj). Using the structures (^) of F'"' and Gf"" into (^) shows that all 
explicit terms in ( |A7|) are zero. Indeed, by multiplying by r^~^ni the term dp{r~^)dQ{r~^) 



in ( [A6a| ) or ( |A6b| ) and projecting on STF tensors, we get a series of terms of the type 
^B-e-2k-2^^ where is a positive or zero integer and i = p + q+ l — 2k. By applying A~^ 
and A~^ we obtain r^^^^'^^fii/ Di 2{B) where the denominators are respectively Di{B) = 
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{B-2i-2k){B-2k + l) and D2{B) = (B -2£-2k + 2){B -2£-2k){B ~2k + l){B -2k + 3). 
When p + q > 1 (which imphes i + k > 1) neither Di{B) nor D2{B) vanish at B = 0, and 
when p + q > 2 (which imphes i + k > 2) D2{B) does not vanish either. Furthermore, by 
multiplying by r^~^ni the second term dp{r~^)dQ{r) in ( |A6 b| ) , projecting on STF tensors 
and applying A~^, we get r^^^^'^^^'^fii/ D^i^B) where the denominator is D^{B) = {B — 21 — 
2k + 2){B — 2k + 3) which does not vanish at 5 = when p + q>2. Since all denominators 
-Di, D2-, -D3 take non-zero values at 5 = we conclude that all terms will be zero at 5 = 
thanks to the explicit factors B present in (|A7|) . Thus, 



can(2) 



0(8,7). (A8) 



The equation ( |A3| ) then shows that for n = 2 the function is 0(8) while 5^, Cl and 
are 0(7). Hence we can write, from 



gl%2) = c^5a(^0(-^))+0(7), (A9a) 



Cn(2) = --Cr ' - Ce^a,^a ['D^i ' ) + ^(8) , (A9b) 

Cn(2) = 0(7) . (A9c) 

Thus it remains to control three terms involving antiderivatives of vectors Ci and and 
having low multipolarities £ = 0, 1. We know that the dependence on of a term with 
multipolarity (i in (3'^^n(2) + + 12 ~ ^) (see e.g. (3.23) in [0), where and £2 cire 

the number of indices on the two moments Vl^ and Vl^ composing the term (notation of 
§ |1V CI). Now one of the two moments is necessarily non-static since for stationary metrics 



the "g-part" of the metric is zero (appendix C in [^), thus li > 2 say. On the other hand, 
to form a vector Ci or Di one needs £2 = thus i.i+i.2 > 3. This, together with the fact 

that £ < 1, shows that the remaining terms in ( |A9|) are 0(7) at least. Thus we have proved 



Cn(2)=0(7,r,7). (AlO) 

We finally deal with the cubic case n = 3. In this case we know that the post-Newtonian 
expansion of the source A^^^^j^g-, starts at 0(6,7,8). [The fact that the spatial components 
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ij of the cubic source are 0(8) instead of the expected 0(6) is not obvious but has been 
proved in [^] — see the proof that A = on p. 424 in |21[]; indeed a possible term 0(6) 
would be made of three mass-type multipoles Ml.] Thus the divergence (^3) with n = 3 



is at least r^an(3) ~ 0(7,8), from which we deduce gcan{3) ~ 0(6,7,8). On the other hand 



the dependence in of a term in g^j^^^g) is 0(8 + S£j — i) (see (3.23) in [Q) and we know 
that S£j > ^ — s by the law of addition of angular momenta, from which we deduce also 
^can(3) ~ ~ = 0(8, 7, 6). These two results imply 

Cn(3) =0(8,7,8). (All) 



Equations ( |A10| ) and (|A11|) are the ones which are used in the text. 



APPENDIX B: THE CONSERVED 2-PN TOTAL MASS 



We first obtain the expression of the total conserved mass at the 2-PN approximation. 
The equation of continuity at this level of approximation reads as 



dt [a (l - 4P/c^)] + dj [aj (l - 4P/c^) 



\idta,j - adtV) - -AaU,d,U + a^udjUk) + 0(5) 



(Bl) 



where our notation can be found in ( p.l|) -(2.4) and (K:.14 )-( 4.15 ). By integrating this equation 
over the three dimensional space we obtain 



d 
di 



1 4 
d-'yia- -ajj - -aP 



J d'^y i^-^adtV - ^^{crU.d.U + a,fe6',f/,)| + 0(5) . 



(B2) 



Several transformations of both sides of this equation yield the equation of conservation of 
mass at the 2-PN approximation, namely 



d 
di 



+4 ( aU' + 2aiUi - 4aj,U - ^dtadtX 



0(5). 



(B3) 
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We now show that the expression ( ^4.21|) we obtained in the text for the general mass- 
type source moment 1^ reduces when £ = to the conserved mass appearing in the square 
brackets of (P3|). When £ = the expression ( [4.21[ ) becomes 



I = YVb=o J Alyl'^l^x + ^{(r,,U - aP) + ^IyI'^'^ 



4 A 



1 + ^ (T,+ 



(dkU[diUk - dkUi] + -dtUdiU 



+ 



120c4 

1 



3 



'P„alU - 2UAd,U + 29.(7j9j(7, - ^(^tU)' - Ud^;U \ + 0(5) . (B4) 



We shall not write down the rather long calculation but only indicate its main steps. One 
must use the equation of continuity ( |B1| ) at the 1-PN approximation, i.e. 



dta + djaj = ^{dtOj, - adtU) + 0(4) 



(B5a) 



together with the corresponding equation of motion 



ai[l + AU/c' 



Oij (1 + AU/c' 



adiV + -[adtUi + ajidjU, - diU,)] + 0(4) . 



(B5b) 



Thanks to the lemma ( ^4.2|) we know that the non-compact supported terms in (p^ ) can 
be freely integrated by parts as if the analytic continuation factor |y|^ and the finite part 
prescription were absent, and as if all surface terms were zero. Let us quote here a list of 
identities, valid up to the required precision and up to the addition of a total divergence, 
which are used in the reduction of 



:PiMU 



1 



120c4 



\y\'dfa 



30c4 
1 



-^\Y?Vid^(Ji — ^|yP<9? 



SttG 



diUd^U 



yidt[dkU{diUk - dkUi)] 



16 
3^ 



Vidt 



1 

4^ 



UdAU 



(B6a) 
(B6b) 
(B6c) 
(B6d) 
(B6e) 
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Using these and other identities, we arrive finally at a manifestly compact-supported expres- 
sion (on which we can remove the analytic continuation factors) which reads as 

+i (<t(7' + 2ffjKi - 4<Tjj(7 - ^dtadtX"j } + 0(5) , (B7) 
in perfect agreement with ( |B3| ) . 



APPENDIX C: COMPUTATION OF THREE CONSTANTS 



To compute the three constants K2, and Kg appearing in ( 4.32 ) one needs to implement 
the construction of the external metric for the interacting multipoles M x Mij (case of ^2), 
M X Mijk (case of ^3) and M x Sij (^2)- The more general cases of interacting multipoles 
M X Milt) and M x Siit) are in fact not more difficult to handle and probably will be 
useful in future work, so we shall compute Ki and n'f for any i > 2. 

In the appendix B of 0, where K2 = 11/12 was obtained, we computed the complete 
M X Mij metric valid all over Dg. Here we shall only compute the terms 1/r (and In r/r) 
in the M x Ml and M x Sl metrics at large distances from the source, since K£ and k,'^ are 
contained in these terms. The quadratic source ( |1.5|) of Einstein's equations, computed with 
the linearized metric ( p.l7|) -( p.l8D and in which we retain only the products of multipoles 
M X M^lt) and M x S^lt), is made of a series of terms of the type dp{r~^)dQ{r~^ F{t — r / c)) , 
where the function F{t) is some time-derivative of a moment M^lt) or SL{t) and where the 
number of space-derivatives acting on is at most two, i.e. p = 0, 1 or 2 (and where q is 
arbitrary). Thus we need only to compute the leading term when r +00 of the (finite 
part of the) retarded integral of dp{r~^)dQ{r^^F(t — t/c)) when p = 0, 1 or 2. When p = 
we know from our previous papers that this retarded integral involves a tail. A computation 
using (2.26) in g and (4.24) in pSf leads to 



X 



2rci JO 

X 



In 



2r 



it---- 

c c 

In r 

2~ 



k=l 



+ 



(CI) 
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(Note that the sum S^^]^(l/fc) is multiphed by a factor 1 in the present formula ( |C1|) and by 
a factor 2 in the formula (2.26) in H].) When p = 1 or 2 there are no tails but the calculation 



is in fact somewhat more complicated. One must use (4.26) in |^3[| to get the polar part at 
i? = of some integrals. The results are 



FP 



B=O^R 



r''d,{r-')dQ[r-'F [t - - 
^ 2(^ + 1) ^'"'"'^''^^^°^"'^''^^ 



(t 



+ 



(C2) 



FPs^oD^' r^a,,(r-^)aQ (r-^F it 



[riij + 5ij)hQ 



2(g+l)(g + 2) 



q+l 



(C3) 



(We denote Q = ai ■ ■ ■ Qg, and in the last term of (|C3|), j means that j has to be excluded from 
the STF operation ( ).) It is then straightforward to obtain the needed nonlinear sources 
with the help of ( |1.5| ) and to apply the formulas (Pl|)-([C^) on each terms of these sources. 
The "p-part" of the metric obtained in this way (see ( |A1| )) is found to be divergenceless up 
to order O (In r/r^) and thus, as a short reasoning shows, the "g-part" of the metric (see 
(|A^)) vanishes at this order. In the case of the interacting multipoles M x Ml, we find the 
metric 



lOO 

'''can(2) 



'''can(2) 



8 ulMM^l^^^ 

e\{i-i) r 

-2 n,^MMi'+^^ 
£\{£ + l) r 
8 riL-iM 



8 ulM 



dxM 



{i+2) 



(C4a) 



+ 



3£ + 4 nL-iMMlr-\ 



e\ e{e + i){£ 



dxM, 



(£+2) 



iL-1 



o 



Inr 



(C4b) 



1,1-3 

"'can(2) 



4 UijLMM^ 



(^+1) 



4 bf + 10£ + 8 nL-mMM^l^'J^ 



+ 



£\{£ + 2) r ' £\ £{£ + l){£ + 2) r 

4 2£^ + 5^ + 4 dj^riLMM^L^^^ 8 2e + h£ + A n^-aMM^+^^a 



l)(£ + 2) 
nL-2M r _,__.^(£+2) 



iiL-2 



£(£+!)(£ + 2) 

£-2 



(C4c) 
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(with c = 1). When C = 2 this metric reduces to (B4) in the appendix B of [Q. In the case 
of the interacting multipoles M x we find 

lnr\ 



1^00 

'<'can(2) 



"'ca.n{2) 



o 



(C5a) 



'^(2) 



i + l)\{£ + l) r 

U UaL-lSiabM 



dxS, 



(£+2) 



+ 1)! r JO 

IQi naL-i{iej)abMSf^L_^ 



bL-l 
(^+1) 



+ 



hir 



(C5b) 



i+l)\{£ + l) r 
lQ{i - 1) naL-2eab{iMSf)^ll^ 



e + i)\{e + i) r 

m naL-2eab(iM 



"'•^'-'j)bL-2 



(C5c) 



In (|C4|) and ( P5|) the moments are evaluated at t — r/c in the instantaneous terms, and at 
t — r /c — x/cin the tail terms. From the metrics ( P^) and ( |C5|) we immediately deduce the 
values of and entering the tail terms in the radiative moments. These are 



k=l ^ 



+ !)(£ + 2) 

£-1 



1 



+ 1) 



fc=i 



(C6) 
(C7) 



We thus find the values quoted in ( |4.33| ). Note that the constants ki and depend on the 
coordinate system which is used, namely the harmonic coordinate system. For instance, the 
constants would he ki + \ and + | in a (perturbed) Schwarzschild coordinate system. 
Note also that 



lim [Ki-\nl] = C 
lim [4-ln£] = C 

f— ► + 00 



(C8a) 
(C8b) 



where C = 0.577. . . is Euler's constant. This can be of interest since the combinations 
K,£ — Ini — C and k,'^ — Ini — C arise in the phase of the Fourier transform of the waveform 



see (3.5) in ||3^ for the mass-quadrupole case i = 2). 
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